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We apply the decoherence formalism to an interacting scalar field theory. In the spirit of the 
decoherence literature, we consider a "system field" and an "environment field" that interact via 
a cubic coupling. We solve for the propagator of the system field, where we include the self- 
energy corrections due to the interaction with the environment field. In this paper, we consider 
an environment in the vacuum state (T = 0). We show that neglecting inaccessible non-Gaussian 
correlators increases the entropy of the system as perceived by the observer. Moreover, we consider 
the effect of a changing mass of the system field in the adiabatic regime, and we find that at late 
times no additional entropy has been generated. 
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I. INTRODUCTION 
A. Outline 

We aim to apply the decoherence formalism to an interacting quantum field theoretical model. The main idea in 
the framework of decoherence ([H-iij], f° r reviews see 041]) is that a macroscopic system cannot be separated from its 
environment. The conventional strategy is to assume the existence of a distinct system, environment and observer. If 
the observer and the environment are weakly coupled, we are allowed to integrate the environment out to study its 
average effect on the system perceived by the observer. Alternatively one could say that the environmental degrees 
of freedom are inaccessible to the observer. This averaging process is an intrinsically non-unitary operation, which 
consequently gives rise to an increase in entropy of the system. A quantum system with a large entropy, in turn, 
corresponds to an effectively classical system. 

It is however difficult to realistically apply the decoherence machinery to quantum field theory: this requires in 
general involved out-of-equilibrium, finite temperature, interacting quantum field theoretical computations. It is of 
course widely appreciated that entropy can be generated as a result of an incomplete knowledge of a system. We thus 
need to keep in mind what quantities are actually measured in quantum field theory: all information in a system that 
can in principle be observed by a "perfect observer" is contained in the n-point correlators of the system. However, 
realistic observers are not capable of measuring irreducible n-point functions of arbitrary order as they are limited by 
the sensitivity of their apparatus. Therefore, it is important to realise that inaccessible higher order correlators, from 
the observer's perspective, yield an increase in entropy of the system. We thus propose the following viewpoint when 
applying the decoherence program to quantum field theory 1 : 

Neglecting observationally inaccessible correlators will give rise to 
an increase in entropy of the system as perceived by the observer. 

As an example, consider some interacting quantum field theory where information is stored in either two-point 
or Gaussian correlators or in higher order, non-Gaussian correlators. The latter are generated generically in any 
interacting field theory. If we assume that the information stored in these non-Gaussian correlators is barely accessible 
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Older work can already be interpreted in a similar spirit [TcL IllH . Here, we propose it as a strict procedure \vX |13|I of how to study 
entropy generation in the 2PI formalism, also see [14| . Recently, during the last stages of writing this paper, we found an interesting 
article by Giraud and Serreau [l5l | addressing the question of entropy production in an interacting field theory from a similar perspective. 
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in experiments, then neglecting this information will give rise to an increase in the entropy. From the Gaussian 
correlators, we can fix the entropy uniquely [Til. [l6j. As before, a quantum system with a considerable amount of 
entropy corresponds to a classical system 2 . We emphasise that this definition can be improved if e.g. three- or 
four-point correlators are accessible through experiments such that knowledge of these correlators is included in the 
definition of the entropy It is important to stress that this procedure does not require a non- unitary process of 
tracing out environmental degrees of freedom. 

In order to apply these rather abstract ideas to a scalar field toy model, let us outline our paper. We will consider 
the following interacting scalar field theory: 

S[<f>, X ] = I d D xC[^x\ = / d D xC Q m +Co[ X ] + Ant[<M , (1) 

where: 

Co[<l>] = -\d^{x)d v 4>{x)^ v - ^ml(t)<t> 2 (x) (2a) 

£o[x] = -^,x(^xWr^™'x 2 W (2b) 

Ant[0,x] = -^(^-hx^x^x), (2c) 

where rj^ = diag(— 1, 1, 1, • • • ) is the D-dimensional Minkowski metric. Here, 4>{x) will play the role of the system, 
interacting with an environment x( x )i where we assume that A > ft such that the environment is in thermal equilibrium 
at temperature T. In this paper, we study an environment of temperature T — 0, i.e.: an environment in its vacuum 
state and we postpone the finite temperature corrections to a future publication. We assume that (4>) = = (%), 
which can be realised by suitably rcnormalising the tadpoles. 

Another application of our calculation is baryogenesis in an early Universe setting where the system is driven out-of- 
cquilibrium by a changing mass term m|(t) generated by a time dependent Higgs-like scalar field during a symmetry 
breaking. For electroweak baryogenesis, we can neglect the Universe's expansion during the phase transition and our 
assumption to work in Minkowski spacetime is well justified. 

As we study out-of-cquilibrium quantum field theory, we work in the Schwingcr-Keldysh formalism. The 2PI (two 
particle irreducible) effective action then captures the effect of perturbative loop corrections to the various propagators 
lA^ and iA x . Of course we will discuss these equations of motion in greater detail in the main text, but when we 
omit all indices and arguments, they have the following structure: 

(d 2 - m 2 )iA^ - J M^A = i5 D (3a) 
(d 2 - m 2 x )iA x - J M x iA x = i5 D , (3b) 

where and M x are the corresponding self-masses. These two equations are non-Gaussian due to the coupling of 
the two fields with coupling constant h. Multiplying equation (|3"a|) by A < p — (d 2 — mj,) _1 (5 D and equation (|3b[) by 
A , x = (d 2 — m x )^ 1 5 D and integrating, one gets the following Schwinger-Dyson equations: 

. + —.'"*•,— 

(4a) 




2 Our definition of classicality differs from the approach used in some of the literature, where for example coherent states with large 
occupation numbers are also considered to have many classical properties, even though these states have vanishing entropy. 
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The dressed ^-propagators (double solid lines) can be written as the sum of the tree level propagators (solid lines) 
and the self- mass corrections due to interaction with the dressed x-field (double dashed lines), and vice versa for the 
environment field. 

Let us at this point explicitly state the two main assumptions of our work. Firstly, we assume that our observer 
is only sensitive to Gaussian correlation functions. This implies that we use only these correlators to calculate the 
entropy 

[3. Secondly, we neglect the backreaction from the system field on the environment field, i.e.: we assume 
that we can neglect the self- mass corrections due to the </)-field in equation (|4b[) . This assumption thus implies that 
the environment remains at temperature T = 0. In particular, the self-mass corrections to the system propagators in 
equation (I4a|) are now given by: 




The first diagram contains the leading order self-mass correction, at order 0(h 2 /ui 2 ), where = rr& + k 2 . It 
effectively is a Gaussian correction to the ^-propagators as it acts just like a known source to 4>. Note that in equation 
([5]) we dropped a diagram at order 0(h 2 X 2 /w^) since it is irrelevant for the argument presented below. 

Looking at equation ([5]) we see that these assumptions are well justified by perturbative arguments provided there 
are no secular effects: the backreaction of the system field on the environment field and the first intrinsically non- 
Gaussian correction occurs only at order C(/i 4 /oj^), which can be appreciated by examining the second Feynman 
diagram on the righthand side of equation ([5]). The 3-point function of the system, an intrinsically non-Gaussian 
correlator, is at one-loop level also perturbatively suppressed at order (D(h 3 /la^). 

Finally, let us just mention that the concept of the pointer basis is frequently discussed in the decoherence literature. 
The pointer basis of our theory in the highly squeezed limit is the field amplitude basis fli| , occurring for example in 
cosmological perturbation theory. Intuitively, this can be appreciated as follows: the Hamiltonian of a squeezed state 
is dominated by the potential term. A system in interaction with an environment at temperature T minimises its free 
energy F = H — TS. The system will realise this by increasing its entropy S, mainly due to increasing the spread in 
momentum (it 2 } since that hardly affects the Hamiltonian, whereas it does significantly affect the entropy. In other 
words: the field amplitude basis is robust under the process of decoherence, qualifying it as a proper pointer basis. 
Note that is a pointer basis only in the statistical sense, such that there is a well defined probability distribution 
function from which a measurement is drawn. 

Having discussed the setup of our theory, the assumptions used and their justification, let us direct our attention 
to discussing potential applications. Our results are relevant for several research areas: the study of decoherence of 
cosmological perturbations, of out-of-equilibrium quantum field theory and of baryogenesis. 



B. Decoherence of Cosmological Perturbations 

Although applications of the framework of decoherence are mainly directed towards experimental efforts (for example 
related to the increasingly relevant field of quantum computing), it also concerns quantum field theory (see e.g.: |17|). 
Research efforts arc primarily focused on addressing the fundamental question of the decoherence of cosmological 
perturbations, see e.g.: [1(1 Ila - l3a |. One of the most important consequences of the inflationary paradigm is that 
it provides us with a causal explanation of how initial density inhomogeneities can be laid out on super Hubble 
scales that seed the large scale structure we observe in the Universe today in for example clusters of galaxies. The 
decoherence formalism applied to cosmological perturbations aims at describing the transition between the quantum 
nature of the initial density inhomogeneities as a consequence of inflation and the classical stochastic behaviour as 
assumed by large scale structure theory. 

In the literature, specific models for example assume that during inflation the UV (or sub-Hubble) modes of 
a field, once integrated out, decohere the IR (or super-Hubble) modes because the former modes are inaccessible 
observationally ( j37H39j , however also see 14011 ) . A similar split of UV and IR modes has been made in the context of 
stochastic inflation, see e.g.: (4fTl46| . In |47| vacuum fluctuations decohere the mean field, turning it into a classical 
stochastic field. In [f6| it is argued that self-mass corrections to the equation of motion for the statistical propagator 
can be rewritten in terms of a stochastic noise term that in turn decoheres the system. In [3l| it was shown that 
isocurvaturc modes decohere the adiabatic mode. 
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C. Non-equilibrium Quantum Field Theory 

In recent years, the study of non-equilibrium quantum field theory has become more and more tractable (for review 
articles see: [3 |4|| ) . A central ingredient in performing these studies is the two particle irreducible action, from 
which quantum corrections to prop agators can be investigated. Out-of-equilibrium X(j> 4 (x) -theory has extensively 
been studied in for example [50l458| . The dynamics of non-equilibrium fermions has been addressed in e.g.: [59| . 

An interesting study has been performed in (60| . where one also studies, under certain assumptions, the dynamics 
of a system field that interacts via a cubic coupling with a thermal bath, which we also consider. Their thermal 
bath consists of two scalar fields with different masses. Very recently, another interesting calculation for A0 4 (x) self- 
interaction has been performed in [l5j where one calculates a decoherence parameter and thermalisation of an initial 
pure state. 

Calzetta and Hu consider in [oil . [(32j also an out-of-equilibrium A<j) 4 (x) theory. What we would refer to as "Gaussian 
von Neumann entropy" is referred to as "correlation entropy" in [62|. They prove an 7?-theorem for a quantum 
mechanical 0(N) model. 

Renormalising the Kadanoff-Baym equations is a subtle business. In A0 4 (x)-theory it has been examined in different 
contexts in j63l47C| . We will also come to address the question of renormalising our cubically interacting field theory. 
Our main finding is that the structure of the renormalised equations of motion differs from the unrenormalised 
equations, which has in general to our knowledge not previously been considered in the literature. 

Furthermore, imposing initial conditions at some finite time to results in additional infinities that have to be 
renormalised separately according to the authors of [7ll - [74l |. Another interesting study has been performed by [lH 
in which the renormalised Kadanoff-Baym equations, again in A</> 4 (a;), are numerically integrated by imposing non- 
Gaussian initial conditions. We differ in our approach as we consider the memory effects from the interacting theory 
at times before io- We can then impose appropriate Gaussian initial conditions at to without encountering initial time 
divergences. 



D. Baryogenesis 

This work is in part inspired by fundamental questions concerning the problem of entropy in field theory, and 
in part by electrowcak scale baryogenesis. The problem is to calculate axial vector currents generated by a CP 
violating advancing phase interface of a true vacuum bubble at the electroweak phase transition. These currents then 
feed in hot sphalerons, thus biasing baryon production. The axial currents are difficult to calculate reliably, since a 
controlled calculation would include non-equilibrium dynamics in a finite temperature plasma in the presence of a 
non-adiabatically changing mass parameter. In this paper we neither include a plasma at finite temperature (this 
will be done in a future companion paper) , nor do we consider scattering of fermions on a non-adiabatically changing 
phase interface. Yet there are important similarities between the problem we address here and baryogenesis: our 
interacting scalar field model © mimics the Yukawa part of the lagrangian of the standard model, whereby one scalar 
field plays the role of the Higgs field, while the other is a heavy fermion (top quark or a chargino of a supcrsymmetric 
theory). The role of the axial current is taken by the entropy which are both sensitive to quantum coherence and the 
phase interface is a time dependent mass parameter rr&(t). The importance of quantum coherence in baryogenesis 
is also treated in (76l - [78j , where a coherent mixture of fermions has been used to generate baryons in grand unified 
theories during preheating after inflation. However, the authors of [76l - [78j treat the interactions phenomenologically 
in the relaxation time approximation. 

Quantum mechanical scattering on bubble walls in a thermal bath may become the dominant mechanism for baryon 
production when the walls are thin, and has been addressed in several papers in the mid 1990s [79H84| mostly in the 
context of baryogenesis within the standard model. Currently, the consensus is that so far no satisfactory solution 
to the problem has been found. Recently Herranen, Kainulainen and Rahkila [85l - [87j have reinvigorated interest in 
the problem, which has gained on timeliness by the upcoming LHC experiments. Their approach is based on the 
observation that the constraint equations for fermions and scalars admit a third particle shell at a vanishing energy, 
fco = 0. The authors show that this third shell can be used to correctly reproduce the Klein paradox both for fermions 
and bosons in a step potential, and hope that their intrinsically off-shell formulation can be used to include interactions 
in a field theoretical settings, for which off-shell physics is essential. The authors have studied both fermionic [HI, |86| 
and bosonic f87| quantum mechanical reflection in the presence of scatterings. However, not all scatterings implied 
by the Kadanoff-Baym equations are taken into account. One important motivation for the present paper is to work 
within a set of approximations where all relevant terms in the Kadanoff-Baym equations arc kept. 
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II. ENTROPY AND PROPAGATORS 



A. The Statistical Propagator and Entropy 



There is a connection between the statistical propagator and the Gaussian entropy of a system also see: 

(Tol [Jj|). In quantum field theory, one can calculate many propagators, with different properties associated with 
each, but not all of them are independent. In this work, we will be primarily interested in solving for the statistical 
propagator of the system. Let us mention that the information contained in the statistical propagator is also encoded 
in the two Wightman functions. Generically, in the presence of quantum fluctuations, one needs complete knowledge 
of the causal propagator in order to solve for the statistical propagator. In the simple free theory example we consider 
in appendix \K\ we can directly solve for the statistical propagator however and no prior knowledge of the causal 
propagator is required. 

The statistical propagator describes how states are populated and is in the Heisenberg picture defined by: 



F<j,(x;x') = -Tr 



= -Tr 
2 



(6) 



given some density matrix operator p(to). The causal propagator roughly describes the number of accessible states 
and is given by the commutator of the two fields: 



iA c 4> (x;x') = Ti p(*o)[0(*),#r') 



Tr 



p(t )MxMx') - 4>{x')^{x)) 



In spatially homogeneous backgrounds, we can Fourier transform e.g. the statistical propagator as follows: 

F^{k : t,t')= [ d(x-x')F^(x;x')e-^< s - s '\ 



(7) 



(8) 



which in the case we will consider in this paper only depends on k = \\k\\. It is only the statistical propagator and its 
various time derivatives that determine the entropy. In short, the entropy is fixed by the area in phase space A the 
state of the system occupies and is given by: 



A 2 k {t) = 4 F(k, t, t')dtdt>F(k, t, t') - {d t F(k, t, t')Y 



(9) 



Throughout the paper, and in particular in this equation we set % = 1. We also set c = 1. The entropy per mode 
then follows as: 



A fc (t) + 1 , {A k (t) + 1 
Sk{t) = log 



A fc (t)-1 {A k {t)-1 
— i — log 



(10) 

Finally, it is interesting to note that the phase space area can be related to an effective phase space particle number 
density per mode or the statistical particle number density per mode as: 



n k {t) 



A fc (t) - 1 



(11) 



In appendix [A] we illustrate our ideas by studying a non-trivial exact case: quantum scattering due to a changing 
mass in the free case, i.e.: the interaction coefficients h and A in equation (|2c|) are switched off. For a free scalar field 
with a smoothly changing mass term, we show that A k (t) = 1 and hence no entropy has been generated by the mass 
change. Secondly, we point out that the reader should not confuse the statistical particle number density in equation 
(fTTjl with the parameter |/3fe| 2 characterising non-adiabaticity of the mass change in equation (|A25bj) . which in the 
literature is often referred to as a particle number as well j88j . This parameter is non-zero, and possibly large, simply 
because the asymptotic in and out vacua differ. 



B. Propagators in the Schwinger-Keldysh Formalism 

The material included in this subsection may well be familiar to the experienced reader, but we include it nevertheless 
for pedagogical reasons and in order to clearly establish our notation. Let us consider the expectation value of an 
operator Q(t) in the Heisenberg picture, given a density matrix operator p(to): 



(Q(t)) = Tr p(t )Q(t) 



Tr 



p(t ) \ Texp 



to 



dt'H(t') } Q(t Q ) < Texp I —i / dt'H(t') 



(12) 
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Figure 1: Schwinger-Keldysh contour with finite initial 
time to and final time tf. 



Figure 2: Schwinger-Keldysh contour where the initial 
and final times in figure [7] have been extended to nega- 
tive and positive infinity, respectively. 



where to < t denotes an initial time, T and T denote the anti-time ordering and time ordering operations, respectively, 
and H (t) denotes the Hamiltonian. If Q in the Scrodingcr picture depends explicitly on time, we should replace Q(to) 
by Q s (t). 

The Schwinger-Keldysh formalism, or closed time path (CTP) formalism, or in-in formalism, is based on the original 
papers by Schwinger [89| and Keldysh [90j and particularly useful for non-equilibrium quantum field theory (also see: 
p)l| - l95| ). According to the CTP formalism, the expectation value above can be calculated from the in-in generating 
functional in the path integral formulation: 

Zi4,jtJ$,J*,p(to)] (13) 
r Vd>+Vc/>oVx^Xo <4,Xo Wo)\<t>o ,Xo) 1^ V<f> + V<t>-S[ct> + (tf) - ^"(*/)] /*° V x +V x ~ 5[ x + (t f ) - 



xexp 



J-D-X, 



where the Lagrangian is given in equation ([T}. We can use the well-known Schwinger-Keldysh contour depicted in 
figures [T] and [2] It runs from to up to tf, where both times can in principle be extended to negative and positive 
infinity, respectively (as depicted in figure [2]). As we will come to discuss, the two paths are not equivalent in an 
interacting quantum field theory, where memory effects play an important role. In this paper, we will extend to to 
negative infinity at some point, but let us for the moment keep it finite. Clearly, these contours are closely related 
to the two evolution operators in equation (|12[) . The fields <j> and X and their corresponding sources and J x 
split up on the upper (+) and lower (-) parts of the contour, where necessarily the conditions <j> + (tf) = <fr~(tf) 
and X + (tf) = X ~(tf) apply. These conditions are indeed enforced by the two functional ^-distributions. The first 
functional integrals in (|13p are over the initial configuration space at to, where the system is specified by the density 
operator (density matrix) p{to)- The path integrals in (|13[) run over the Schwinger-Keldysh contour in figure [1] or [21 
Expectation values of n-point functions are obtained by varying the generating functional (|13[) as follows: 



Tr 



p{to)T[<j>{ Xl ) . . . <t>{x n )]T[j>( yi ) . . . <j>(y k )} 



S n + k Z{J lP (t )] 



i5Jt{xi) ■ ■ ■ i5JZ(x n )i5Jl(yx) ■ ■ ■ i5Jl(y k ) 



provided that < tf and < tf for all j, and where J = ( J^, J^). We can now define the following propagators: 



iA++(x;x') 

iA-+(x;x') 
iA+-(x;x') 



= Tr 



Tr 



= Tr 



Tr 



p{to)T[4>{x')4>{x) 
p(t )T[4>(x')$(x) 

p(t )4>(x)4>(x') 

p(t )$(x')$(x) 



= Tr [p(to)<f> + (x')<f> + (x) 
= Tr [p(to)$-(x')$-(x) 
p(t )4>-(x)$ + (x') 
p(to)r{x')4>+{x) 



S 2 Z[J,p(t 



i5J*(x)i5jf(x') 
S 2 Z[J lP (t )] 



= Tr 



Tr 



i8r{x)i5Jl{x') 



S 2 Z[J,p(t )} 



j=o 



,7=0 



i5JZ{x)i5Jl(x') 
S 2 Z[J,p(t )] 



i&jUx)i5Jt{x' 



,7=0 



,7=0 



(14a) 
(14b) 
(14c) 
(14d) 
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We define the various propagators for the x _ field analogously. In the absence of a condensate for x an mixed two 
point functions, such as: 

(n\4>( x ')x(x)\n) , (15) 

vanish by virtue of the interaction term (|2c|) . In equation (fT4|) . iA^ + (x; x') = iA^(x;x') denotes the Feynman or 
time ordered propagator and iA~^~(x; x') represents the anti-time ordered propagator. The two Wightman functions 
are given by iA^ + (x; x') and iA^~ (x; x 1 ). From equation (P) we infer that the free Feynman propagator obeys: 

V x iA++(x; x') = {d 2 x - m 2 )iA++{x- x') = i5°(x - x') , (16) 

where d x = rj^d^d^, and where the same identity holds for zA^q(x; x'). In the presence of interactions the equation 
of motion for the Feynman propagator becomes much more involved and we will discus it shortly. One can easily 
show that at tree level the Wightman functions obey the homogeneous equation: 

V x iA+~{x- x')=0 = V x iA~+{x- x 1 ) . (17) 

Identical relations hold for the free x _ P ro P a g a tors. The four propagators defined above are not independent. The 
Wightman functions for example constitute the time ordered and anti-time ordered propagators: 

iA++{x-x') = 6{t-t')iA^ + (x;x') + 6(t' -t)iA+-(x;x') (18a) 
iA^-(x;x') = 6(t' -t)iA^+(x;x') + 6(t-t')iA+-(x;x') (18b) 

where t = x°, t' = x a> , and where this identity holds for the x _ P r opagators as well. Appreciate that: 

iA+ + (x;x') + iA^-{x;x') = iA^+{x; x') + iA+-(x; x') (18c) 

iA-+(x;x') = iA+-(x';x), (18d) 

are exact identities and they are also satisfied by the x-prop a g a tors. We write the four Green's functions in the 2x2 
Kcldysh propagator matrix form: 

d9) 

which at tree level obeys, 

V x iQ^{x- x 1 ) = ia 3 S D (x - x') , (20) 
where a 3 = diag(l, —1) is the third Pauli matrix, which we can also write as: 

(a 3 ) ab = aS ab , (21) 

where a, b = ±. 

Let us define some more Green's functions. In subsection III Al we already defined the causal and statistical prop- 
agator, but let us for completeness list them again. The causal Green's function, also known as the Pauli- Jordan or 
Schwingcr or spectral two-point function, iA c ^ = tA^ J = A$ = p^, is given by: 

lAfa; x') = Tr (p(t ) <£(*')]) = iA~ + (x; x') - iA+~{x- x') , (22) 
and the statistical or Hadamard two-point function, = A^ , is given by: 

F^x;x>) = ±Tr[p(t ){t(x'),$(x)}] = \ (iA^ + {x; x') + iA+~{x; x')) . (23) 

The retarded («A r ) and advanced (?A a ) propagators are defined as: 

iA^(x;x') = tA^^x^-iA+'^x') 

= -[iA--{x ] x')-iA~ + {x ] x')]=9{t-t')iA%{x ] x') (24a) 

iA%{x;x') = iA++(x;x')-iA^+(x;x') 

= -[iA--(x;x')-iA+-(x;x')] = -6(t' - t)iA%{x; x') . (24b) 
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Moreover, we can express all propagators lA 1 ^ solely in terms of the causal and statistical propagators: 

iA+- (x; x') = F^x; x') - \iA%{x; x') (25a) 

iA-+(x;x') = F^x-x')+ l -iA%{x-x') (25b) 

iA++(x;x') = F^x;x') + ^sgn(t-t')iAl(x;x') (25c) 

iA--(x-x') = F^x^-hg^t-tyAl^x'). (25d) 

Since F^ — F^ and (lA^y — — iA^, the relations above correspond to splitting the various Green's functions into 
their hermitian and anti-hermitian parts (for that reason we do not put an i in front of F$). The definitions of the 
retarded, advanced, causal and the statistical propagators and the relations between them easily extend to the x-field. 

C. The Kadanoff-Baym Equations 

In order to study the effect of perturbative loop corrections on classical expectation values, one often considers the 
effective action. In this subsection we will calculate the 2PI effective action, using the Schwinger-Keldysh formalism 
outlined above. The 2PI effective action is the relevant functional to consider because it captures the interaction 
of the <f>- and x-fields in the right way. Varying the 2PI effective action with respect to the propagators yields the 
so-called Kadanoff-Baym equations that govern their dynamics. These equations of motion contain the non-local 
scalar self-energy corrections or self-mass corrections to the propagator. 

In the present subsection, we shall mainly follow [48], Hj, [9a, [9?| ■ We can extract the Feynman rules from the 
interaction part of the tree level action (|2"c)) : 



A„t[0,x] = -$> Q(X») 3 + ^(x»)V(*)) . (26) 

a— ± ^ * ' 

The Feynman propagator is promoted to tQ^ and each vertex has two polarities: plus (+) and minus (— ), such that 
the minus vertex gains an extra minus sign +i\ as compared to the standard perturbation theory plus vertex — iX. 

The 2PI effective action can be obtained as a double Legendre transform from the generating functional W for 
connected Green's functions with respect to the linear source J and another quadratic source (see e.g.: [48j|). In the 
absence of field condensates the background fields vanish: 

4> a = (n\4> a \n) = o (27a) 

X a = (n\x a \ty = 0, (27b) 

in which case the variation with respect to the linear or quadratic sources can easily be related. In particular, the 
definitions of the four propagators in equation (TH| remain valid. The effective action formally reads [48l. [93l [98l [991] : 

r^rVA^A^] = S^^^ + lTrln^Afr^+'-TrHizAf)- 1 ] (28) 



l T . ^ 2 <S[0 a ,X a ] A ab , l T 5 2 S[(j) a ,X a ] .ab , r (2)r7a -a a ah Aai 



Here, r^ 2 ^ denotes the 2PI contribution to the effective action. Moreover, we omitted the dependence on all variables 
for notational convenience. Several Feynman diagrams contribute to the effective action of which the one and two 
loop order contributions arc given in figure [3l We can now write the effective action up to two loops as: 

T[iAf,iAf] = T [iAf ,iAf] +T 1 [iAf ,iAf] +T 2 [iAf ,iAf] (29) 
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Figure 3: Contributions to the 2PI effective action up to two loop order. The double solid 
lines denote (^-propagators, whereas the double dashed lines correspond to x-P ro P a 9 a t° rs - 



where the subscript denotes the number of loops and where: 

T [iAf,iA x b ] = ( d D xd D x' ^(dl-ml)6 (x-x')S ab iAf(x';x) (30a) 

a,b=± 

+ f d D xd D x' ^(dl-rn 2 x )6 D (x-x')6 ab iA b x a {x';x) 

a,b=± 

T^Af^} = ~Trln[%A^(x;x)]-^Trhi[t^(x-,x)] (30b) 
T 2 [iAf,iA ab ] = — d D xd D x' Y, ab(iA x b (x';x))+— d D xd D x' £ ab (iA ab {x;x')) iAf{x';x),(Wc) 

a.b—i: a.b— ± 

where Tr denotes a trace over both spacctimc variables and the Keldysh indices ±. The equations of motion for the 
propagators result as usual from the variational principle: 

1 -4 ' . x 1 = (31a) 

\abl 

h — = 0. (31b) 
Explicitly, they yield: 

l^-ml^-^-I^Af^xOl-V^^ 6 ^;^) 2 - (32a) 
a -{&l - m x )8 D (x - x')5 ab - I [iA x b (x- x')] ~ l + ^ab (iA ab { X] x')f + l ^-abiA ab {x; x')iAf(x; x') = . (32b) 

We will bring these equations into a more familiar form by multiplying by 2atA b £(x' ; x") and 2aiA b x {x' ; x"), re- 
spectively, and then integrating over x' and summing over b = ±. This results in the following one-loop Kadanoff- 
Baym |X OOj ] equations for the elements of the Keldysh propagator iQ(x;x'): 



6T[iAf,iA x b ] 



{d 2 x - ml)iAf(x; x 1 ) ~J2 C d ° x i M !p c (^ Xi)iAf{ Xl \ x') = aS ab i5 D (x - x') (33a) 
(pi - m x )iA ab {x-x') - Y c i & D xiM* c {x\ xi)iA cb (x x ;x') = a5 ab i5 D (x - x') . (33b) 



c=± 

E 

c=± J 

where the self-masses at one loop have the form: 

iMftex!) = ~2ac ^J^y = -V (^ c (*;*i)) 2 (34a) 
iM™(x; Xl ) = -2ac J^-±J—£l = -— (zA^x; Xl )f - ih 2 iA a x c (x; x^iA^x; x x ) , (34b) 

where in the last step we used the hermiticity symmetry of the operator iQ, according to which, iA ac (x;x') = 
iA ca (x' ; x). The Feynman diagrams contributing to the one- loop self-mass are given in figure 0] We have chosen the 
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Figure 4'- Contributions to the self-masses up to one-loop order. Again, the double solid lines 
denote <j> -propagators, whereas the double dashed lines correspond to x-P ro P a 9 a ^ ors - Hence, the 
first two Feynman diagrams contribute to the self-mass of x( x )i an ^ on hj the third diagram 
contributes to the self-mass of 4>(x) . 



definition of (|34|) such that the structure of the self-mass resembles that of the propagators. The factor 1/2 in (|34|) 
originates from the symmetry factor of the one-loop self-mass diagram. 
Equation (|33a[) consists of the following four equations: 



(^ 2 -m^A+ + (.T;.T') 
(^-m|)*A+-(z;z') 
(d 2 x -ml) l A- + (x;x') 
{dl-ml)iA--{x;x') 



iM+~ 

4> 


'(x;y)iA 4 


- + {y;x')_ 


= i8 D {x-x') 


iMf- 

4> 


'{x\y)iA~ 4 


-(y;x')~ 


= 




'(x;y)iA 4 


- + (y;x')_ 


= 


iM7~ 
4> 


'{x\y)iA 4 


-{y;x') 


= -i5 D (x-x 



(35a) 
(35b) 
(35c) 



but in the light of equation (|18[) . only two of them are independent. Note that we have another set of four equations 
of motion for the x- field. In the end, we will be interested in solving this equation of motion in Fourier space, e.g.: 



iAf(x;x>) 



W) 1 



Tl Af(k,t,t')e 



/\ ik(x-x') 



iAf{k,t,t') = / ^-'{x - x ! )iA a 4 b {x-x')e- l ^ s - s '^ 



Such that equation (|35|) transforms into: 



{d 2 t +k 2 +ml)iA++(k,t,t')+ / dti iM++(k,t,t 1 )iA+ + (k,t 1 ,t') - iM+-(k,t,h)iA^ + (k,h,t') 



(d 2 + k 2 +ml)iA+-(k,t 7 t')+ / dti iM+ + (k,t,t 1 )iA+-{k 7 t l7 t')-iM+-(k,t 7 t 1 )iA--{k,t 1 ,t') 



(36a) 
(36b) 



(37a) 
iS(t - 1') 
(37b) 



{d 2 + k 2 + m 2 )iA 4 + (k,t,t')+ / dti iM- + (k,t,t 1 )iA+ + (k,t 1 ,t')~iM 4 -(k,t,t')iA- + (k,t 1 ,t') 



— OO 

oo 







(d 2 + k 2 + m 2 )iA 4 -(k,t,t')+ / dti iM 4 + (k,t,t 1 )iA+-(k,t 1 ,t')-iM 4 -(k,t,t 1 )iA 4 :~(k,t 1 ,t') 



(37c) 

(37d) 
-i8(t-t'). 



Note that we have extended to — > — oo in the equation above. Again, we have an analogous set of equations of motion 
for the x-field. In principle we can solve these coupled equations of motion only numerically in full generality. Our 
strategy is to push the analytical calculation forward as far as possible, before relying on numerical methods. Before 
we make an important simplifying assumption, let us first consider the renormalisation of our theory. 



11 



III. RENORMALISING THE KADANOFF-BAYM EQUATIONS 



In order to renormalise equation of motion (|35p or (|37[) above, we need to Fourier transform also with respect to 
the difference of the time variables: 

{ c\ D k 

ab(w>\ = / _ r/\«b( h n^ik-(x-x) 



iAf(x;x') = J _*A?(fc'')e**- 
iAf(k^) = (& D (x-x')iAf(x]x')e- lk ^ x - x '\ (38b) 



There is a subtlety: for the moment we neglect the time dependence in the mass term. We only use this assumption 
to renormalise. In the end it turns out that we need a mass independent counterterm to cancel all divergences in 
our theory, which allows us to consider a time varying mass term again. In fact, as we assume there is no residual 
dependence on the average time coordinate (t + 1')/2 in iA^ fc (fc M ), equation (f38|) coincides with a Wigner transform. 
Fourier transforming equation of motion (|35p yields: 

(-k^-ml- l M+ + (k^A+ + (k^+ i M+-(k^ + (k^ = i (39a) 

{-k^-ml-iM+ + ^))iA+~{k^)+iM+-{k^A--{k^) = (39b) 

{-W-ml + iMj- (kn>^ + (k^-tM^ + (k^)tA+ + (k^ = (39c) 

(-k^-ml + iM^-(k»))iA^-(k»)-iM^+(k»)iA+-(k») = -i, (39d) 

Here and henceforth, we use the notation fc M fc M = — fcg + ^ 2 to distinguish the four- vector length from the spatial 
three- vector length k = \\k\\. Because of the convolution, the equations of motion above are local in Fourier space. 
Let us remind the reader again that analogous equations hold for the x~P r0 P a g a t° rs - As already announced in the 
introduction, we shall not solve the dynamical equations for both <fi- and x-propagators. Instead, we shall assume the 
following hierarchy of couplings: 

h < A (40) 

and expand the solution in powers of h/X <C 1. In fact, we shall solve the system only at order (/i/A)°. This does 
not imply that the hy 2 ^ interaction is unimportant: we will only assume that A is large such that the %-field is 
thermalised by its strong self-interaction. This allows us to approximate the solutio ns of the dynamical equations for 
X as thermal propagators which we derived in the appendix in equation (|B10[) . see jl0l| : 



»A++(fcM) = — — „ +27rS(k^ + m 2 x )n^(\k \) (41a) 



\ 



lA ^ {kn = k^ + m l + l e +2nS{k ^ + < )nC " ilk ° l) (41b) 
iA+-{k») = 2iTS(k^ +m 2 x )[0(-k°) + n^{\k°\)] (41c) 
iA-+(k») = 2TrS{k^ + m 2 x )[6{k°)+n^(\k°\)] , (41d) 

where the Bose-Einstein distribution is given by: 

"?(*") = ^ = ^- (42) 

with fcs denoting the Stcfan-Boltzmann constant and T the temperature. Let us remark that assumption (|40[) allows 
us to compute the quantum corrections to the (^-propagators as it depends solely on x-propagators running in the 
loop. We neglect the backreaction of the system field on the environment field, such that the latter remains in 
thermal equilibrium at temperature T. This assumption is perturbatively well justified as already discussed in the 
introduction. Furthermore, we neglected for simplicity the C(A 2 ) correction to the propagators above that will slightly 
change the equilibrium of the environment field. Note finally that, in our approximation scheme, the dynamics of the 
system-propagators is effectively influenced only by the usual 1PI self-mass correction. 

In this paper, we consider only an environment field x m its vacuum state at T = and we postpone the finite 
temperature corrections to a future publication. Any divergences, if present, originate from the vacuum contributions 



12 



to the self- masses, i.e.: the vacuum propagators at T = are a useful case to consider anyway: 

iA++{k») = — — 9 (43a) 



k^k^ + to 2 , 
l 



IL 



A x"(* M ) = ttt^t-t: (43b) 



k^k^ + to 2 , + ie 



iA+-(k") = 2n5{k ll k f " + m 2 x )9(~k°) (43c) 
iA-+(k>*) = 2Tr6(k fl k fl + m 2 x )6(k°). (43d) 

We evaluate the Feynman self-mass iM^ + (x; x') following from equation (|34aj) where we make the simplifying as- 
sumption m x — > 0. Let us briefly expatiate justifying this assumption as at first sight it seems that m x — > makes 
our approximation scheme more susceptible to undesired backreaction effects 3 . It is a priori not at all clear that the 
backreaction is negligible: if we examine the second Feynman diagram on the right-hand side of equation (j4b| we see 



that the leading order backreaction occurs at order 0(h 2 /w 2 ,). Since in our setup u> x = k 2 + m 2 , = k 2 , it is clear that 
the backreaction on deep IR (infrared) Fourier modes of the environment field is pcrturbatively unsuppressed. Despite 
that, it does not spoil the perturbative arguments employed in the introduction: the influence of the environment field 
on the system field is still perturbatively under control. In order to see this, let us consider the first non-Gaussian 
contribution on the right-hand side of equation ((SJ). Indeed, one can show that the IR part of the inner loop is phase 
space suppressed: the IR part of this integral is given by A 4 q[(k v k u — mV)(q a + fc CT )(q ,<T + k a )}~ 1 ~ fl 2 /m| when 
nitf, > h ~ p. Nevertheless, we admit it would be worthwhile to examine these integrals for m x ^ and see whether 
the results presented in this paper are robust under this change. 
We thus need to evaluate: 

l Aj+(x-; x') = J jgL iA J+(fc^(— ') . (44) 

This integral can be performed in arbitrary dimensions by making use of two straightforward contour integrations 
and [l02l . 1l03j : 

d"- 1 /? r - If 00 Jd-3 (kx) 

:x f(k) = —^r Mk D - \> „_ 3 /(fc), (45) 



-e 



(27T) - 1 ( 47 r)^ Jo [\kx) 

which is valid for any function f(k) that depends solely on k = \\k\\. J^(kx) is a Bessel function of the first kind. This 
yields: 

i A++(x; x') = 1 2 > J . (46) 

* 47nr Ax ++ {x;x') 

Here, Ax^_ + (x, x') is one of the distance functions between two spacctime points x and x' frequently used in the 
Schwingcr-Kcldysh formalism and given by: 

Ax% + (x,x') = -(|t-i'|-^) 2 + ll^-^'H 2 (47a) 

Ax 2 + _(x,x') = -{t-t' +ief + \\x -x' || 2 (47b) 

Ax 2 _ + (x,x') = - (t-t' -ief + \\x-x" || 2 (47c) 

Ax 2 __{x,x') = -{\t-t'\+ief + \\x-x' || 2 . (47d) 

We thus immediately find from equations Q34ap and (f4"6)) : 

,1,2 r2 ( B _n i 

iM++{x-x') = - — -^ ~^r^ ■ (48) 

4> K ' J 2 16it d Ax 2 + D + - 4 {x;x') V ; 



3 We thank Julicn Serrcau for this useful comment. 
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The other self-masses can be obtained from this expression using the appropriate e pole prescription as indicated 
in equation (|47p . We will now rewrite this expression slightly in order to extract the divergence. For an arbitrary 
exponent /3 ^ D, f3 7^ 2, we can easily derive: 

1 1 „ 9 1 



d 2 



Ax%Ax;x>) 08 - 2)08 -D) Ax^ + 2 (x; x>) 



Furthermore, recall 104L 1 10 



(49) 



1 _ 47T 2 ^ 

A<; 2 (* ;a; ') ~r(^) lx xj ' 



Let us also recall the similar identities for the other distance functions: 



d 2 
d 2 
d 2 



1 



Ax D _- 


l 2 (x;x') 




1 




l 2 {x-x>) 




1 



4-7T 2 

W) 



id D (x-x') 



Ax°- + 2 (x;x') 

We now arrange equation (j48|) . using (|49"]) and (|50a|) : 




0. 



(50a) 

(50b) 
(50c) 
(50d) 



2-r2( D 



iM++(x;x') = - 



iti 2 T 2 ( 



1 



64tt d 



(D-3)(D-4) 



9* 



47T 2 



^ D + -'(i;i') A< + 2 (x;x') iW) 



-p = 2vi(5 r> (x — a;') 



(51) 

Here, the scale fi has been introduced on dimensional grounds. If we Taylor expand the term in curly brackets 4 around 
D = 4, we find: 



iM++(x;x') 



ih 2 T{§ - 1)m d - 4 
16ttt(.D-3)0D-4) 



i5 D (x-x') 



ih 2 



1287T 4 



r<) 2 



\og{li 2 Ax 2 {x-x')) 



Ax 2 ++ (x\ x') 



+ 0{D-4) 



(52) 



We have been able to separate a local (D — 4) 1 divergence and a non-local finite term to the self-mass. In order to 
precisely cancel the divergence, we can thus add a local counterterm, i.e.: an ordinary mass term of the form: 



iM±t t (x;x')= T 



ih 2 T(§ - 1) M D - 4 
167T* (L>-3)(L>-4) 



i5 n (x-x'). 



(53) 



The relative sign difference of iM^ ct (x; x') is due to equation (|50h>|) . We are left with the following renormalised 
self-mass: 



lh 2 T 2 { D_ x) 



64tt d (D-3)(D-4) 



d 2 



1 



2n- B 



(x;x') Ax°- + 2 (x;x>) 



We will now perform a spatial Fourier transform in order to solve for the dynamics this term generates: 



^ + r + en(M,0 



\x-x')iM++ en {x;x')e 



(54) 



(55) 



By introducing a regulator in order to dispose of the overall surface terms arising from two partial integrations, we 
can easily convert the partial derivatives. Using several analytic extensions, we obtain: 



D-13 _ -D + l 



k^{D -3)(D-4) 



3-0 1 , D 



7T2 fc 2 



r(^) 

k D -\i\At\ +e)^ 
2 D - 4 r(D-3) 



-jfc(|At|-ie) 



(56) 



-if D -5 (fc(i|At| +e)) 



4 Note that in the minimal subtraction scheme, one would also expand the term multiplying the Dirac delta function around D = 4, 
which gives rise, once integrated at the level of the equation of motion, to a finite local contribution to the mass of cf>. 
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where k = \\k\\, At = t — t' and where K v (z) is the modified Bessel function of the second kind. We expand this result 
around D = 4: 



*M++ n (fc,i,t') 



__ e - 4 M|A t |-.)( 7E + log 



2tn 2 (\At\ - le) 



(57) 



■y/%\M\ + eduK v {ik{\At\ - it)) 



v=l/2 



0{D - 4) 



Here, 7e is the Euler-Mascheroni constant. Moreover, the scale n introduced earlier combines nicely with the other 
terms to make the argument of the logarithm dimensionless as it should. Indeed, we need to find an expression for 
the derivative with respect to the order v of K v (z). Starting from the general expansion: 



fc=0 



1 



T(k-v+l)k\ \2 



% \ 2k— v 



l 



Y(k + v + \)k\ \2 



Z \ 2k+v 



we immediately derive: 

d v K v (z) 



v=l/2 



e 2 [Chi(2z) - Shi(2z)] 



7e + log(2z) 



(58) 



(59) 



dt- 



cosh t — 1 



t 



dt- 



sinhi 



t 



where Chi(2z) and Shi(2z) are the hyperbolic cosine and hyperbolic sine integral functions, respectively, defined by 
the expressions on the second line. In our case, the variable z is imaginary, so it proves useful to extract an i and 
convert this expression to the somewhat more familiar sine and cosine integral functions, defined by: 



si(z) 
ci(z) 



dt 
dt 



smt 
cost 



We finally arrive at: 



64fc7T 



-'ifc|Ai| 



7E + log 



+e 



2in 2 (\M\ - it) 
(ci(2fc(|Ai| - it)) - iai(2k(\M\ - it)) 



(60a) 
(60b) 

(61) 



+ 0(D-4), 



where we have set the e regulators in the exponents to zero as the expression is well defined. Rather than going several 
times through the calculation above to determine the other self-masses, we make use of a few analytic extensions. 
Observe for example that if At > 0, Ax ++ (x, x') and Ax- + (x,x') coincide, hence the expressions for self-masses 
iM^ cn (k,t,t') and iM^ + (k,t 7 t') should also coincide in that region. All we need to do is to sensibly analytically 
extend to At < 0. We will thus need: 



si(— z) = — si(z) — 7r . 

If A< < 0, we have to carefully make use of the t pole prescription in the cosine integral function: 



ci(-2fc(-A£ + it)) 



dt- 



cost 



-2/c(-At+ie) 



t 



dt 



-2fe(-At) 





poo 


/ dt + 


/ dt 


J —ML 


J2fe(-At) 



cost 



(62) 



(63) 



- log(ze) + log(-te) + ci(2fc(-Ai)) 
-j7r + ci(2fc(-Af)) . 



We thus find the following expressions for the renormaliscd self-masses: 

iM^ cn (k, t, t') = (d 2 t + k 2 ) iZf{k, t, f) . 



(64a) 
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whore: 

Z^(k,t,t') 

Z^(k,t,t') 



64fc7T 2 

h 2 

64kn 2 



e^ kAt ( 7 e + log 



2^ 2 \At 
k 



7T 

T i— 
^ 2 



,±ifc|At 



(^ci(2fc|At|) =f zsi(2fc|Ai|)^) 



(64b) 



2^ 2 |Ai| 



T «|sgn(Ai)J + e ±?feA * (ci(2fc|Ai|) T zsgn(Ai)si(2fc| A«|) 



Firstly, appreciate that zM^ + (fc, i, t') and iM^ (k, t, t') need not be renormaliscd. The reason is that these expressions 
do not contain a divergence in D = 4, which can be seen from equations (|50c| - |50d|) . Moreover, the local counterterm 
which we add to renormalise iM^~ [k 1 1, t') contains the opposite sign as compared to iM^~ + (fc, t, t') because of equation 
(I50bl) , which we already stated in equation (|53| . Finally, we have sent all e regulators to zero as the expression above 
is well defined in the limit At — > 0. 

We performed two independent checks of the calculation above. Firstly, one can renormalise via a calculation in 
Fourier space (rather than position space) . We show that the two results agree in appendix [C] Secondly, one can 
calculate the retarded self-mass directly from the position space result using (|52l) and compare with the result obtained 
from We show that the two results agree in appendix ID] 

If one were to evaluate the two time derivatives in the expressions above, one would find a divergent answer in the 
limit when At —> 0. We also show this in appendix[C] This does not reflect an incorrect renormalisation procedure. It 
is crucial to extract the two time derivatives as presented above in order to properly take the effect of the self-masses 
into account as only now Z^ b (k,t,t') is finite at coincidence At — s- 0. Indeed, this is most easily seen in position 
space 5 . 

Let us compare these expressions with existing literature. In e.g. flU [66| it is derived that the renormalised 
equations for A</> 4 (a;) theory have an identical structure as the unrenormalised equations. In our theory, clearly, the 
structure of the two equations changes as we need to extract an operator of the form (d 2 + fc 2 ), as derived in equation 

(El. 



IV. DECOUPLING THE KADANOFF-BAYM EQUATIONS 

Having renormalised our theory, we are ready to massage the Kadanoff-Baym equations ([57)) in two different 
ways. Firstly, we will write Kadanoff-Baym equations in terms of the causal and statistical propagator such that 
they decouple. This is of course a vital step required to solve the Kadanoff-Baym equations in the next section. 
Secondly, we show that when we write the equations in terms of the advanced and retarded propagators, the one-loop 
contributions preserve causality as they should. 

Note that the structure of the self-mass (IM1) is such that we can construct relations analogous to equation (|25l) . 
which of course hold identically for \: 

Z+-(k,t,t') = Zl{k,t,t')-^Zl{k,t,t') (65a) 

Z^ + (k,t,t') = Z$(k,t,t') + ±iZ%(k,t,t') (65b) 

Z++(k,t,t') = Z^{k,t,t') + ^gn{t-t')iZl{k,t,t') (65c) 

Z^-(k,t,t') = ZF(k,t,t')-l S gn(t-t')iZ%(k,t,t'), (65d) 



5 One can easily recognise that the structure of the renormalised self-masses in equation II64H is identical to the d'Alembertia n in Fourier 
space. The presence of Z^ b (k, t, t') induces time dependence in the propagator. A similar phenomenon has been observed in |106(I . where 
this phenomenon is referred to as a "finite wave function renormalisation" , in which the effect of gravitons on fermions in an expanding 
Universe is investigated. 
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such that we find from ([64]) : 



l r 
2 



z7+(fc,t,t')+^r(fc,i,t') 



64fc7T 2 



cos(fcAt) 7e + log 



2/* 2 |At| 



• ci(2fc|Ai|) ) +sin(fc|At|) (si(2/c|Ai|) - | 



z+-(k,t,t')-z- + (k,t,t') 

h 2 



(66a) 



(66b) 



64fc7T 2 



2 cos(fcAi)sgn(Ai) (si(2fc|Ai|) + |) + 2sin(fcAi) ( ci(2fc|At|) - 7e - log 



" k ' 


)] 


.V|At|_ 





The expressions for x differ due to We can derive a system of two closed equations for the causal and statistical 
propagator by adding and subtracting equations (|37cp and (|37b[) . In order to obtain the equation of motion for the 
causal propagator (JT]), we subtract (|37bl) from (|37cl) to find: 



(% + k 2 +ml)iA c Jk,t,t') + - 



dh 



[iM^+(k, t,h)- iM+~{k, t, ti)} B gn(ti - t') 



tA c Jk,h,t') = 0. (67) 



Using equations (|64|) and (|65j) we find: 

(d 2 + fc 2 +m 2 )A^(fc,M') - (a 2 + k 2 ) J f dt 1 Z^k,t,t 1 )A^(k,t u lf) = 0. (68) 

Note that equation (|68[) is causal, in the sense that no knowledge in the future of the maximum of t,t' is needed to 
specify iA^(k, t, t'). Moreover, at one loop the evolution of «A^ requires only knowledge of the Green's functions in 
the time interval between t' and t, and is thus independent of the initial conditions at to = — oo. Finally note that 
we deleted the i in front of iA^ in the equation of motion above to stress that A^ is real to prepare this equation for 
numerical integration. 

In order to get an equation for the statistical Hadamard function ([6]), we add equation (|35b[) to (|35c[) to get: 



(d 2 + k 2 + m 2 )F (M,O 
1 



(69) 



dh 
dh 



iM 



r++ 



(k, t, h) - iM+-(k, t, h) + iMT+{k, t, h) - iM7~{k, t, h) F^k, h,t') 



0,renv "i "J ""^ V v l "> "IJ ' if, V™) "> " ivJ 0,renV 

iM+~(k, t, h) + iM-+(k, t, h)} sgn(h - t') - iM++ n (k, t, h) - iM^- e Jk, t, h)] iA%{k, h,t') = 0. 
Again using (|64l) and (I65|) we find the relevant differential equation for the statistical propagator: 



{d 2 + k 2 + m l)F ttl {k,t 1 t l )-{d 2 + k 2 ) 



f dhZ^hMF+fah,?) - f dhZ^{k,t,h)A%{k,h,t') 

J — oo J — oo 



0. (70) 



We can thus say that the equations of motion for the causal and statistical propagator have decoupled in the following 
sense: the differential equations haven been brought in triangular form. Note that equations (|6"5|) and ([70|) together 
with the causal and statistical self- masses in equation (|66[) . represent a closed causal system of equations suitable 
for integration in terms of an initial value problem. Given the knowledge of F and iA c for both \ anc l 4>i &U other 
Green's functions can be reconstructed from equation (|25p . This strategy was used (see [48[ and references therein) 
to study the dynamics of out-of-equilibrium quantum statistical (scalar and fermionic) field theories. Indeed, we will 
solve equations (|6"5|) and ([70)1 numerically in the next section. We emphasise however that the form of equations ([70)1 
and (|68[) differs from the ones found in [48j . The renormalised equations of motion have a different structure than the 
unrenormaliscd ones, which is not taken into account in e.g. [48l |66|. 

Before doing so, let us show that the one-loop self-masses do not spoil causality in another way: the retarded and 
advanced Green's functions only receive information from the past and future light cone, respectively. Now subtracting 
equation (|37b[) from (|37a|) one obtains: 



(d 2 + k 2 + m 2 )iA%{k, t, t') + / dhiM^ ren (k, t, h)iA\{k, h,t') = i5{t - t') . 



(71) 
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Making use of equation (|D4[) we find: 

iM; iIen (k, t,h) = - {d 2 + fc 2 ) 0(i - h)Z%% t, h) 
Equation of motion (|71j) transforms into: 



(a 2 + k 2 + ml)i^(k, t, o - (a 2 + fc 2 ) / dtiz^fc, t, ti)*A r (fc, *!, o = - t') 



(72) 



(73) 



The retarded self- mass gets contributions only from within the past light cone, i.e.: when h < t. 

Similar to equation (|73[) . we can subtract equation (|37c[) from (|37al) to obtain the equation of motion for the 
advanced propagator: 



(a 2 + fc 2 + m 2 >A^(fc,M') 



dfi 



(»A^+„(fc, t,ti)- iM^+(k, t, h)) iA+ + (k 7 h,t') 
- (iM+-(k,t,h) - iM^- ea (k,t,h)) iA-+(k,h,t') 



This yields: 



(a 2 + fc 2 + ml)iA%(k, t, t 1 ) - (a 2 + fc 2 ) / dhz;(k, t, hYA%{k, t u t') = iS(t - 1') 



iS(t - t') . (74) 



(75) 



where we find an analogous relation for the advanced self-mass: 

iMl Ien (k,t,h) = - (d 2 + fc 2 ) 9(h - t)Z%(k,t,h) 



(76) 



As expected, iM^ Ien (k,t,h) acquires contributions from the future only, i.e.: when h > t. Rather than solving for 
the causal propagator, we could alternatively solve for the retarded propagator or the advanced propagator. We will 
however not pursue this in the present work. 



V. NUMERICALLY SOLVING THE KADANOFF-BAYM EQUATIONS 



Let us once more explicitly write down the equations of motion of the causal and statistical propagators (|68|) and 
([70]) we will numerically tackle in this section: 



(a 2 + fc 2 +m 2 )A£(fc,M')- (a 2 + fc 2 ) / &hZl{k,t,h)&%{k,h,t') =0 (77a) 



(a 2 + fc 2 +m 2 )i^(fc,M')-(3 2 + fc 2 ) 



dhZ;(k,t,h)F4k,h,t') - / dhZ^(k,t,h)^l(k,h,t') 



= 0. (77b) 



The causal and statistical self-masses are given in equation (|66p . In particular, we will be interested in two cases: 



i^4>(t) = m a = const 

m 2 Jt) = A + B tanh(pt - 50) , 



(78a) 
(78b) 



where we let A and B take different values. Let us take a closer look at the two equations of motion above. Clearly, 
we first need to determine the causal propagator. Note that equations of motion (|77j) depend on two variables, i.e.: 
for each t', we have to solve this equation of motion 6 . The self- mass corrections contribute only through a "memory 



Alternatively, we could have written down the equations of motion of the causal and statistical propagator where the operator acts on 
the other leg of the propagator, on t' . Then we would have to solve these four equations of motion simultaneously. Needless to say the 
two methods are completely equivalent. 
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kernel" (memory integral over time) between t' and t. The boundary conditions for determining the causal propagator 
are as follows: 



A C (M) = 
d t A%{t,t')\ t=v = -1 



(79a) 
(79b) 



Condition (|79a[) has to be satisfied by definition and condition (|79b[) follows from the Wronskian normalisation 
condition due to the commutation relations. 

Once we have solved for the causal propagator, we turn our attention to the second equation (|77b[) . Suppose we 
would not have sent to — > — oo. The equation for the statistical propagator then would have been of the following 
form: 



{d 2 t 



ml)Ft{k,t,t')- (d 2 + k 2 



f dhZ^t^F^tut') - f to x z$(k,t,hWt{k,hJ) 

to J to 



(80) 



Clearly, equations (|80|l and (|77bp are not equivalent. Equation (|77bp contains a memory kernel from the infinite 
past up to t and t', whereas equation (|80|) only contains memory kernels from to onwards. This corresponds to an 
interaction that is switched on non-adiabatically at time to- To understand this, consider replacing the coupling 
constant h hidden in the self- masses (l66l) with 7 : 



h -> h6(ti - t ) • 



(81) 



The step-function would then have transformed equation (|77b[) to (|80[) which mimics switching on the interaction 
between the two fields at some finite time to- The two standard Schwinger-Kcldysh contours presented in figures 
[T] and [5] are thus not equivalent in interacting quantum field theories where memory effects play an important role. 
Alternatively, we could say that non-locality, generic for any interacting quantum field theory, enforces the memory 
kernel to start at the in finit e past. This effect has, in the context of electromagnetic radiation, been recognised and 
investigated by Serreau [l07l | . In the work of Borsanyi and Reinosa (67], HH the memory integral, extended to negative 
infinity, plays an important role too. They suggest to use that in connection with a generalised dissipation-fluctuation 
theorem. 

Needless to say, we have to start at some finite time in our numerical analysis. We therefore make the assumption 
to approximate the propagators in the memory kernels from the negative past to to with the free propagators: 



(d 2 t +k 2 +ml)F (t) {k 1 t,t')-{d 2 + k 2 ) 



dt 1 Z%(k,t,t 1 )F* ee (k,t 1 ,t') + / dtiZ£(M,ii)F (Mi,O 



'(1 



d^(fc,Mi)A;< frco (Mi,0- f dt 1 Z| , (A,* > ti)AJ(fc J ti > t / 

J to 



(82) 



=0, 



where Ff*>(k, h, t') and Af lee (k,tx,t') are the free propagators obtained in equation (|A9[) . This approximation 
induces an error of the order 0(/i 4 /u;4). An alternative approach has been outlined in [75j where, for X<f> 4 (x) theory, 
non-Gaussian initial conditions at to are imposed. We can explicitly evaluate the infinite past memory kernel: 



M| 0G 



(k,t,t',t ) = (d 2 + k 2 ) 
h 2 



t 



dti 



Z%{k, t, h)F^(k, h,t') - Z$(k, t, h)Af ee (k, h,t') 



32ir 2 uj m J _ 
We change variables to r = t — ti to find 



*° cos(kt + Lu in t' - (k + u- m )ti) 
dti 

DO 



t-tl 



Mp e (k,t,t',t ) 



dr 



cos((fc + Lu in )r - uj in (t - t')) 



(83) 



(84) 



[cos(w in (i - t'))ci({k + uj in ){t - t )) + sm(cj in (t - t'))si((k + uj in )(t - to))] 



7 Note that t > t\ by construction. 
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We postpone the discussion of imposing boundary conditions for i 7 ^ at to to subsection IV Al Equation (|77p transforms 
into: 



{d 2 + k 2 +ml)A%{k,t,t')-(d 2 + k 2 ) J d^(fc,Mi)A^(Mi,i') = 
(d 2 + k 2 + m%)F^k % t, t') - Mp c (k, t, t', t ) 



(85a) 
(85b) 



(d 2 + k 2 ) 



dt 1 Z^(k,t,t 1 )F (j> (k,t 1 ,t') 



dtiZ| , (M,*i)A$(Mi,* / ) 



We have now prepared the problem for numerical integration. In the numerical code, we take to = and we let 
pt and pt' run between and 100. In order to solve differential equation (|85bj) . we thus need to evaluate two more 
memory kernels. The one involving the causal propagator can be computed immediately Once we have solved for 
the statistical propagator, our life becomes much easier as we can immediately find the phase space area via relation 
(|9]). The phase space area fixes the entropy. 

Differential equation (|85b[) merits another remark. If we let t — > to, we encounter a logarithmic divergence in 
Mjre e (fc, ^ ^ ^ ag (-.j^ ^ log(a;) as x — > 0. This divergence is only apparent. Intuitively, this should of course be the 
case as we introduced the boundary time to by hand and no divergences should arise consequently If h = const, the 
time to is introduced as a fictitious time, hence observables cannot depend on to. Of course, neglecting the memory 
integral from negative past infinity to to introduces a dependence on t . Thus, removing the distant memory integrals 
completely is equivalent to setting h — s- h9(t\—to) as in equation (|81|) . We will prove that this logarithmic divergence is 
only apparent rigorously by rewriting equation (|85a[) for the causal propagator and (|85b>|) for the statistical propagator 
in a different form, and by using the symmetry properties of the propagators. Focussing first on the equation of motion 
for the causal propagator (|85aj) . note that we can transfer the t derivative to a ti derivative by using the fact that the 
causal self-mass (|66b[) is a function of At = t — ti only: 



d 2 J ^ dti^(A I t,t 1 )AJ(fc,ti,0 = -d t J dtiS tl Z$(M,*i)A$(Mi,*' 

= d t j ' dti^(M,*iRA£(Mi/) 



(86) 



= -z;(k,t,f) + J dtiZKM.tO^A^Mi,*'), 

where we partially integrated in the second line (the boundary terms vanish by virtue of equation (|79aj) and Z^(k, t, t) = 
0), and we used the commutation relations in the third. We transform the equation of motion of the statistical 
propagator (|85b|) analogously to find: 



(d 2 + k 2 + r4)A^(k,t,t') + Z;(k,t,t')- dtyZl{k,tM){dl+k 2 )A%{kM,t') =0 (87a) 



ml)Ft(k, t, t') - Mp c {k, t, t', t ) - \d t Z%(k, t, to)F+(k, t , t') - d t Z% (k, t, t )A£(fc, t , t') 
-Z% (fc, t, t') + Z%(k, t, to)dtMk, to, t') - Zl (k, t, to)d to A%{k, to, t') 

dt 1 Z c 4 ,(k,t,t 1 )(dl+k 2 )F4k,t 1 ,t / )+ f dt 1 Z^{k,t,t l )(d 2 tl +k 2 )A%{k,t 1 ,t l ) -0. (87b) 



These two differential equations should be completely equivalent to equation J85J). In fact, we have found a non-trivial 
test of our numerical code: the results of equation (jB"5j) and of the two equations above should agree. We will show 
this in due course. 

Now, we can see another logarithmic divergence appearing in dtZ^{k, t, to) in equation (|87b() when we send t — > to. 
The reader can easily verify that the logarithmic divergences in Mp cc (k, t, t' , to) and dtZ^(k, t, to)F<p(k, to, t') in equa- 
tion (|87b[) above cancel to leave a finite result when t — > to if we set 2i 7 ^(fc,t ,t') = cos(cj in (t — t'))/u- m . We thus 
find that at order 0(h 2 /u^) no divergences at to remain and we expect that a similar treatment would cure these 
type of apparent divergences at higher order: clearly, to has been introduced by hand so this should not lead to any 
irregularities. 

Let us finally make some remarks about the literature. The authors of [7l| - [73 | study out-of-equilibrium \<fi 4 (x). They 
encounter, after renormalisation, a residual divergence in their theory at the surface of initial boundary conditions 
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Figure 5: Statistical propagator in Fourier space for 
a small coupling h/p = 3/2 (black) and a larger one 
h/p — 4 (dashed). Because of a non-zero coupling, 
we observe that the S-function, present in the original 
dispersion relation, has broadened to a "quasi-particle 
peak", roughly of a Breit-Wigner form. If the coupling 
increases, the "quasi particle peak" broadens further. 
Clearly, when h ^> ui^ in the strongly coupled regime, 
we have a "collection of quasi particles". We used 
k/p = 1, m^/ p = 1 and p/p = 1. 



Figure 6: Statistical particle number density rife as a 
function of k/p. We used: h/p = 1 and m^/p — 1 
(black), h/p = 2 andm$/p = 1 (black dashed), h/p = 1 
and m^/p — 0.5 (blue) and h/p — 2 and m$/ p = 0.5 
(blue dashed). In the UV nk vanishes irrespective of 
the value of h/ p or m<p/p. Particles are only produced 
by the interaction in the statistical sense in the IR. 
The mass only influences the IR behaviour. Moreover, 
one can show by appropriate rescalings of the functions 
above that is given by equation in the UV. 



at to which they choose to renormalise separately. We differ in their approach as we do not find these residual 
divergences. The infinite past memory kernel preci sely takes care of these as can be appreciated from the previous 
discussion. This is also the case in the approach of [67], |H[ mentioned before. 



A. Constant Mass Solutions 



The constant mass case is interesting as we can make non-trivial statements based on some analytical calculations. 
The bottom line is that the generated entropy is constant. The argument is rather simple. When to^ = const, we 
have: 

F <p {k,t,t')=F <p (k,t-t'). (88) 
Using a few Fourier transforms (|38[) . we have: 

F*{k,0) = J -^F<p( k n (89a) 

r°° Ah 

Wfc,At)| At=0 = -i J —k°F^) (89b) 
f°° dk a 

d t ,d t Ft(k,At)\ At=0 = / —k 2 MW), (89c) 

The right-hand sides do no longer contain any time dependence. Hence, the left-hand sides are also time independent. 
This implies that the phase space area is constant, and so is the generated entropy. As we insert some finite t , 
we expect to observe some transient dependence of the entropy on time because we approximated the propagators in 
the infinite past memory kernel with free propagators. When this behaviour has died out, the entropy should settle 
to its constant value derived from equations (|89[) . (|9]) and (fTO)) . Indeed, this constant entropy does not necessarily 
equal 0. We interpret this non-zero entropy as the entropy generated by the coupling to the second field, which in 
the effective action acts as a source for F^. Effectively, this opens up phase space for the system field that previously 
was inaccessible for it. More accessible phase space for the system field in turn, implies that less information about 
the system field is accessible to us and hence we observe an increase in entropy. 

In order to evaluate the integrals above, we have derived the statistical propagator in Fourier space in appendix lEl 
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The result is: 



--sgn(fc>(fc 2 



k 2 ) 



k^kv + m 2 . 



h 2 

32tt 2 



Ioe 



2 7E - |||sgn(fc")0(fc o 2 - 



/t 2 

327T 3 



log 



27E)+C s g n W(fco 2 -fc 2 ) 



(90) 



To gain some intuitive understanding, we depicted the statistical propagator for h/p = 3/2 and for h/p = 4 in figure 
[5j When h/p = 0, we have a <5-function dispersion relation as usual but in the presence of a non-zero coupling, the 
5- function broadens to a so-called "quasi particle peak" of a Brcit-Wigncr form. For h/p = 4, this peak is still well 
pronounced, but when we enter the strongly coupled regime, this simple picture breaks down when the resonance 
becomes broad and we can no longer sensibly talk about a "quasi particle" , but rather we should think of a "collection 
of quasi particles" . 

The integrals in equation (|89p can now be evaluated numerically to yield the appropriate initial conditions. For 
example when k/p = 1, m^/p = 1 and h/p = 4, we find: 



F^,(k/p=l,At)\ Ai=0 = 0.35196 
d t F 4> (k/p=l,At)\ At=0 = 
d t >d t F 4> (k/p=l,At)\ At=0 = 0.73120 ; 



(91a) 
(91b) 
(91c) 



Clearly, equation (|91b|) holds for all values of k, m$ and h as the integrand is an odd function of k°, which can be 
appreciated from equations (|89bp and . The numerical value of the phase space area in this case follows from (|9"Tj) 
and © as: 



A ms = 1.01461 > 1 . 



(92) 



Hence also: 



0.04327 > . 



(93) 



where the subscript ms is an abbreviation for "mixed state" . 

From the phase space area A ms we can easily obtain the statistical particle number density (|11[) . It is interesting to 
study its behaviour as a function of k. Figure [6] clearly shows that in the deep UV (ultraviolet) the particle number 
density vanishes: the interaction between the two fields only produces particles (in the statistical sense) in the IR. 
Moreover, using figure El we can show: 



n k {h,m^) 



' 



(94) 



in the deep UV. In fact, we can estimate the constant of proportionality £ appearing in equation (f9"4"|) as C — 0.0008 
which turns out to be insensitive to the value of the mass of the system field and the coupling h. The mass only 
influences the IR behaviour, as expected, which can also be appreciated from figure [5] Note finally that the formal 
divergence of derived quantities, such as the total particle number per volume N/V = J d 3 k / (2tt) 3 n k or the total 
entropy per volume S/V = f d 3 k/(2ir) 3 Sk does not pose any problems for the dynamics we are about to solve since 
these quantities do not enter the equations of motion. 

So far, wc postponed the discussion of imposing boundary conditions for numerically determining the statistical 
propagator. We just proved that, independently of how one imposes initial conditions, the phase space area should 
settle to a constant value and for a specific choice of parameters, we have been able to calculate this constant in 
equation (|92p . One could think of at least two separate ways of imposing boundary conditions: "pure state initial 
conditions" and what we will henceforth refer to as "mixed state initial conditions" . If we constrain the statistical 
propagator to occupy the minimal allowed phase space area initially, we set: 



F<f,(to, t ) 



1 



d t F4t,t )\t=t = 



d t ,d t F^t')\ t=tl=t0 = ^ 



(95a) 
(95b) 
(95c) 
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and Wi n is determined from equation (jAllI) . This yields: 

A k (t ) = l, (96) 

such that Sk(to) = 0. Physically, this means that despite the fact that interactions enlarge the accessible phase space 
of the system field, we force it to occupy a minimal area initially and let it evolve 8 . Alternatively, we can impose 
mixed state initial conditions, i.e.: the values calculated from equations ([89]) and ((901) : 



i^(Mo,*o) = / ( 97a ) 

d t F4fc,*,t')| t=t , =t0 = (97b) 

dfc° 



d tl d t F 4 {k,t,t% =v=to = J_ —klF^), (97c) 
such that initially: 

A fc (to) = A ms = const, (98) 



and also S ms > 0. Clearly, boundary conditions (|97[) can only be evaluated numerically for each choice of parameters. 
Needless to say we are completely free to impose any other type of initial conditions as well, but we consider the two 
cases above to be physically well motivated if the system is close to its minimum energy state. 

In figure [71 we show the phase space evolution for both pure state initial conditions (black line) and mixed state 
initial conditions (red line). For pure state initial conditions, the evolution is precisely as anticipated. The phase 
space area increases from its minimal area Afc(fo) = 1, to the asymptotic value A ms calculated in equation (|92p . For 
mixed state initial conditions, we firstly observe some transient behaviour which eventually decays. We then smoothly 
evolve to A ms . The initial transient is due to our assumption to approximate the propagators in the memory kernel 
from past infinity to to with free propagators. As is apparent from ([84]) . its effect becomes less important as time 
elapses. 

From the evolution of the phase space area (for pure state initial conditions), we can immediately find the time 
evolution of the entropy Sk{t) in figure This shows that entropy has been generated by interaction with an 
environment that is in the vacuum state assuming that some observer is only sensitive to Gaussian correlators. The 
entropy eventually settles to its asymptotic value S ms calculated from A ms and equation (|10p . As can be anticipated, 
the generated entropy per mode is small: both system and environment are in a state close to the minimum energy 
state (T = 0). 

We conclude that when the mass m$ is a constant, no further entropy is generated if we start with mixed state 
initial conditions. However, we do observe a generation of entropy if we start with pure state initial conditions. This 
increase in entropy can be understood by the system's tendency to evolve towards the vacuum state of the interacting 
theory. 

Let us return once more to figure [7J The fact that our numerical asymptote is located slightly above the one 
calculated from (|92p can be attributed to the accuracy of the implementation of the infinite past memory kernel. This 
can be appreciated from figure [5] where we test the accuracy of our code. Clearly, the numerically found asymptote 
decreases towards A ms as accuracy increases. Moreover, observe that the initial violent oscillations in Afc(i) decrease 
as accuracy improves. Also, we have chosen wAi, where At is the step size of the numerical integration, such that we 
resolve all the oscillations. For a N = 2000 run at tp = 100, we have uiAt ~ 0.071 for the parameters used in figure 
[7] Also, we can observe a "beating" phenomenon that persists even if the accuracy increases (and that can hence 
not be attributed to numerical artifacts) . It is caused by a frequency mismatch by approximating the propagators in 
the infinite past memory kernel by free propagators. Finally, let us discuss figure 1101 Here, we show the difference 
between the phase space evolution calculated from equation ([55]) and from ([57)) . The dashed line is more accurate 
than the dotted one, which in turn is more accurate than the solid line. Clearly, the differences disappear when the 
accuracy improves. This confirms our numerical analysis in a non-trivial way. 



If we would not include the infinite past memory kernel Mp ca (k,t,t' ,to) and indeed consider a coupling between two fields that is 
switched on non-adiabatically at some finite time to as previously discussed, the pure state initial condition would be the natural choice 
for this problem. 
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Figure 7: Phase space evolution for constant for 
pure state initial conditions Ak(to) = 1 (black line) and 
mixed state initial conditions (red line) Ak(to) = A ms . 
In both cases the phase space area settles to the constant 
value A ms , indicated by the dashed line and calculated 
in equation S92j) . We use k/p = 1, m^/p = 1, h/p = 4 
and N = 2000. 

A 




Figure 8: Convergence for the phase space evolution 
presented in figure [7j The black and red lines are 
identical to the ones in figure an d most accurate 
(N = 2000,). The gray and yellow lines are calculated 
with N = 1000. The other parameters are kept fixed. 
Clearly, the difference between A ms and the numerical 
asymptotes decreases as the accuracy increases. 



Figure 9: Entropy generation for the system field <f> 
through interaction with the environment \ in the vac- 
uum state. As both fields are in a vacuum state, the en- 
tropy generation is relatively small. We used pure state 
initial conditions Sk(to) = and k/p = 1, m^j p = 1, 
h/p — 4 and N = 2000. The entropy settles to a con- 
stant value Sm S calculated from the value of A ms of 
figure^ and equation flOV . 
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Figure 10: Test of our numerical code. This plot shows 
the difference of the phase space area calculated from 
equation \85\) and \87ty for different values of N but 
the same values for the other parameters k/p = 1, 
vn^/p = 1 and h/p = 4 for pure state initial condi- 
tions. We used: N — 1000 at t = 100 (solid line), 
N = 1000 at t = 50 (dotted line) and N = 1000 at 
t = 25 (dashed line). The difference between the two 
methods disappears as the accuracy of the numerical 
evolution increases. 



B. The Decoherence Time Scale 



We define the decoherence time scale to be the characteristic time it takes for the phase space area Afc(i) to settle 
to its constant mixed state value A ms . One can suppose that such a process is described by a differential equation of 
the following form: 

SA k (t)+T k (h,LJ <f> )SA k (t)=0, (99) 

where SAk(t) = A k {t) — A ms and where T k (h,ui^) is the decoherence rate. This equation is equivalent to h k — 
—T k (n k — n ms ), where n k is defined in equation (1111) and n ms is the stationary n corresponding to A ms . We anticipate 
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Figure 11: Exponential approach to A ms . We study, 
for different initial conditions, differences of Ak(t) on 
a logarithmic scale defined by equation HlOlp . Clearly, 
the decoherence rate does not depend on the initial con- 
ditions. We used k/p — 1, m^/p — 1, h/p = 4 and 



Figure 12: Decoherence time scale. We confirm the 
scaling relation for the decoherence time scale antici- 
pated in equation flOOV . For all plots we took k/p — 1 
and N = 1000, and furthermore we used m^/p — 1, 
h/p = 4 (solid red line), m^/p = 1, h/p = 1 (dashed 
red line), m^/p — 4, h/p = 4 (solid blue line) and 
m 4>/p = 4, h/p = 1 (dashed blue line). 



N = 2000. 



that the decoherence rate depends both on the coupling constant and on the energy of our system field 9 . The 
following intuitive picture is helpful: the solution of equation (|99[) results in an exponential decay to the mixed state 
value 5Ak(t) cx exp[— Tk(h, <rf$)t]. Furthermore, a stronger coupling h should result in a larger value of Tk{h,ui,f,). 
However, a larger energy w'i = mi + k 2 should be reflected in a smaller value of Tk(h, uj^). On dimensional grounds, 
we thus anticipate: 



where 7 = const. Let us now test this expected scaling relation. Looking back at figure [7j we see in the first few 
time steps that Afc(t) oscillates. Clearly, the time scale of these fluctuations has nothing to do with the decoherence 
time scale, but rather can be fully attributed to numerical accuracy. To capture the decoherence time scale correctly, 
we thus consider the difference A™ s (t) — A ps (£) of the evolution of the phase space area A™ s (i) using mixed state 
initial conditions and A ps (i) using pure state initial conditions. On a logarithmic scale, we observe in figure [Til an 
exponential decay towards A ms (solid line). Moreover, the slope should not depend on the particular choice of initial 
conditions. To this end, we also calculate the difference of A" ls (i) — A)J ps (£), where A)J ps (i) follows from setting 
1 < A.^ pb (t ) < A™ s (io) initially. In order to do this, we kept the value of F0(fc,< o ,^o) identical to the value it had 
for the mixed state boundary conditions but reduced the value of dfdtF^k, t, t')\ t _ t ,_ t such that the inequality 
1 < AjJP^to) < A™ s (t ) is satisfied. The resulting decoherence rates precisely coincide as seen in figure [TT1 where we 
plot: 



and likewise for A£ ps (i) (dashed line). 

We can repeat the steps outlined above for a different choice of parameters m and h. By rescaling the obtained 
decoherence rates by a factor of u^/h 2 , we can test the scaling relation (|M| . All decoherence rates now precisely 
overlap as we depicted in figure 1121 We can thus estimate the value of the constant of proportionality 7 appearing in 
equation (|100[) : 




(100) 




(101) 



h 2 

T k (h, U4) = (0.0101 ± 0.0003)— . 



(102) 



Ideally, we should of course take the 0(h? /u)^) to oj^ through the dispersion relation into account. 
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This relation gives the decoherence rate for our particular model. This result is nothing but the single particle decay 
rate: 

lm(tM++) 1 h 2 . 
= , = ^77 ' (103) 

where we have used equation (|C2[) and e.g. |60lll08j . Let us compare the result (|102[) to the literature. Let us remark 
that most of the calculations have been performed in an expanding Universe setting, or with a different model, so 
it is hard to compare this result quantitatively. In [39j it was found that, for a different model during inflation, the 
decoherence rate is proportional to the spatial volume, which we do not find. 



C. Changing Mass Solutions 

Finally, let us discuss the evolution of the phase space area when rr&(t) is changing according to equation (|78b[) . 
The analytic expression for the statistical propagator in Fourier space we previously derived in equation (|90[) is no 
longer valid. Introducing a time dependent mass m|(t), generated by a time dependent Higgs-like scalar field, breaks 
the time translation invariance of the problem. Consequently, the statistical propagator F(k,t,t') no longer depends 
only on the time difference of its time variables At = t — t', because considering rn\{t) introduces a proper time 
dependence on the average time coordinate r = (t + t')/2 in the problem. When the mass is changing rapidly, we 
can only rely on numerical methods. However, asymptotically, where the mass settles again to a constant value, the 
analysis performed in the previous subsection should still apply. 

We impose mixed state boundary conditions as in equation ([97]) such that Afc(i ) = A ms initially. Of course, the 
value of the mass inserted to calculate these initial conditions, is the value of the initial mass, valid before the mass 
jump. 

If the mass changes non-adiabatically, this results in a significant particle creation according to the discussion in 
appendix [A] We can thus identify the following regimes: 

|/?fc| 2 <C 1 adiabatic regime (104a) 

\Pk\ 2 ^ 1 non — adiabatic regime , (104b) 

where j3k is given is equation (|A25b[) . For the parameters we used in figures [T31 and [T5l we are in the adiabatic regime 
(|/3 fc | 2 =3.5-10- 4 ). 

In figure [13] we study a mass increase from m^j p = 0.75 to m^/ p = 2 (black line). This decreases the phase space 
area and consequently the entropy decreases, which we depicted in figure HU Intuitively, a larger mass of the 0-ficld 
reduces the effect of the quantum corrections of the x-field. Hence Ams < Ami, where Ams is the constant phase 
space area calculated for m^/p = 2.0 from three Fourier integrals as in equation (j9"Tj) . Likewise, Ams corresponds to 
the phase space area calculated for m^/p = 0.75. The relevant behaviour to compare with is the constant mass phase 
space evolution for / p = 2 which we also depicted in figure [13] in gray. Clearly, the late time asymptotes of the 
two functions coincide and we conclude that, also at late times, no entropy has been generated. As we are in the deep 
adiabatic regime, this is to be expected. 

Now let us study the opposite: a mass decrease from m^/p = 2 to m^/p = 0.75. This is depicted (black line) in 

figure 1151 Clearly, the phase space area increases from Ams to Ams and we plotted the resulting entropy increase in 
figure [TBI Again, we compared this evolution with the phase space area calculated for a constant mass m^/p = 0.75 
(gray line). As the two asymptotes also coincide in this case, we conclude that no entropy has been generated at 
late times by the mass change. Of course, it would be very interesting to see what happens when we study the same 
process in the non-adiabatic regime and we hope to address this question in a future publication. 

If we compare the evolution of the entropy in time in these two cases with the free case Sk(t) = 0, the interacting 
case reveals much more interesting behaviour. Firstly, due to the presence of an environment field, the constant value 
to which the entropy settles asymptotically is different from zero, unlike the free case. Secondly, a changing mass 
induces dynamics: the entropy depends on time and evolves from one value S'ms to another S^s or vice versa. 



VI. CONCLUSION 



We apply the decoherence framework to quantum field theory. We consider two scalar fields, one corresponding 
to a "system field" , the second to an "environment field" , in interaction via a cubic coupling. Here, we consider an 
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Figure 13: Phase space area decrease due to a mass 
increase from m^/p = 0.75 to vn^jp = 2 (solid black 
line). We used h/p = 4, k/p = 1 and N = 2000. We 
thus observe a slight entropy decrease. The solid gray 
line denotes the constant mass phase space evolution 
for m^/p = 2. As the two asymptotes coincide, we 
conclude no entropy has been generated at late times 
by the mass change. Ami and A„] are the constant 
mixed phase space areas calculated for m^/p — 0.75 
and m<f,/p = 2, respectively. 
s 

0.06 r 




Figure 14: Entropy decrease for the case presented in 
fiaure rnft Because of the mass increase, the phase space 
area decreases which results consequently in a drop in 
the entropy. 




Figure 15: Phase space area increase due to a mass 
decrease from m^/p — 2 to m^j p = 0.75 (solid black 
line). The other parameters are the same as in fig- 
ure \13\ The solid gray line denotes the constant mass 
phase space evolution for m$/ ' p = 0.75. Again the two 
asymptotes coincide and no entropy has been generated 
at late times by the mass change. 




Figure 16: Entropy increase for the case presented in 
figure [751 Clearly, a mass decrease gives rise to a slight 
entropy increase. 



environment in its vacuum state (T = 0) and postpone finite temperature contributions to a future publication. We 
neglect the backreaction of the environment field on the system field, assuming that the former remains at T = 0. 

We advocate the following point of view regarding a sensible application of decoherence to quantum field theory: 
for some observer inaccessible higher order correlators give rise, once neglected, to an increase in entropy Sk of the 
system. This is inspired by realising that correlators are measured in quantum field theories and that higher order 
irreducible n-point functions arc usually perturbativcly suppressed. In this work, we assume that our observer is 
only sensitive to Gaussian correlators and will hence neglect higher order, non-Gaussian correlators. Neglecting the 
information stored in these higher order correlators, gives rise to an increase of the entropy of the system. If the system 
initially occupies the minimal area in phase space (characterised by a pure state with Sk(to) — 0), we numerically 
calculate the evolution of the entropy Sk(t) in figure [Sj Also, we calculate the asymptotic value of the phase space 
area A ms and the entropy 5 ms as a function of the coupling h, the mass and k to which these functions evolve. 
This increase in entropy can be understood from the system's tendency to evolve towards the vacuum state of the 
interacting theory. Even though we do not solve the full 2PI equations at one loop, we have strong numerical evidence 
that within our approximation scheme, the system evolves towards its correct stationary interacting vacuum state. 
If the system, however, starts out initially occupying the state characterised by S ms , we observe no further increase 
in the entropy. Furthermore, we calculate the decoherence rate in equation (|102p which charaterises the exponential 
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rate at which the system approaches its stationary state. 

Secondly, we study the effect of a time dependent mass m^(i) of the system field. Now, calculating Sk(t) can only be 
addressed numerically. Starting out at mixed state initial conditions Sk(to) = S ms > 0, we observe an entropy increase 
(decrease) due to a mass decrease (increase) as depicted in figures [TB] and 1141 respectively. By comparing with the 
constant mass evolution for the entropy, we conclude that no additional entropy has been generated asymptotically 
by the mass jump. As we study mass changes in the deep adiabatic regime, it remains to be investigated whether this 
statement also holds in the non-adiabatic regime. 

We also would like to draw a few somewhat more technical conclusions. It is important to stress that in interacting 
field theories where memory effects play a crucial role, one cannot just insert initial conditions at some arbitrary finite 
time to j because one then neglects the memory effects existing from the infinite asymptotic past to t . In numerical 
computations however, one has to start at some finite time. We therefore approximate the propagators in the memory 
integral from the past infinity to to by free propagators, inducing a perturbatively suppressed error. 

Also, it has not been previously appreciated in the literature that renormalising the vacuum contribution in the 
Kadanoff-Baym equations can actually significantly change the structure of these equations. In order to properly take 
account of the renormalised self-mass contribution we had to extract two time derivatives which can readily be seen 
from equation (|64p . 
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Appendix A: Quantum Effects of a Changing Mass: Free Case 



It is interesting to compare our results to a non-trivial exact case: scattering of a changing mass field in the spirit 
of Birre l l and Davies [H[ . The solutions presented here stem f rom the cosmological particle creation literature (based 
on |l09l . Illd |) and are originally due to Bernard and Duncan 111 | . Let us consider the action of a free scalar field: 

S[<f>] = J d 4 x l—d^^xyrT - \mftt)<j?{x)\ , (Al) 

where as usual ?^ t „ = diag(— 1, 1, 1, 1) is the Minkowski metric, and where we consider the following behaviour of the 
mass m,p(t) of the scalar field: 

ml(t) = (A + Btanh(pt)) . (A2) 

From (jAT|) . it follows that: 

(9 t 2 -^ + m2(t))0(x)=O, (A3) 
The vacuum causal and statistical propagators follow from (O and ((6]) where p(to) = |0)(0| as: 

i&%{x;x') = (0\[fa),$tf)]\0) (A4a) 
F^x-x') - i(O|{0(x'),<^)}|O>. (A4b) 
Let us quantise our fields in D-dimensions by making use of creation and annihilation operators: 

fa) = / (^7 (%Mt>** + ■ (A5) 

The annihilation operator acts as usual on the vacuum: 

%|0) = 0, (A6a) 

and we impose the following commutation relations: 

[ag.ajp] = (27r) I, - 1 (5 D - 1 (fc-A r '). (A6b) 
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Hence the equation of motion for the mode functions <j>k(t) of <j){x), defined by relation (|A5[) . follows straightforwardly 



as: 



(d? + k 2 + ml(t))4> k {t)=0, (A7) 
where k = \\k\\. The mode functions determine the causal and statistical propagators from (|A4p completely: 



i&%{k,t,t') = sin(w(£ - 1')) (A9a) 



*AJ(fc > *,0 = (/> k (t)cj>l{t') - <t>k(t'Wk(t) (A8a) 

F*(fc,M') = i{M*0#(*) + &(*)dE(*0}- (A8b) 

Let us at this point for completeness calculate the constant mass causal and statistical propagators in Fourier space. 
We just insert a constant mass m^, rather than a changing one as in equation (| A2|) : 

i 

F^,{k,t,t') = — cos(o;(t - 0) > ( A9b ) 

where ui 2 = m\ + k 2 . Let us now return to the changing mass case. The physical picture is clear: we would like to 
study reflection and transmission, i.e.: scattering, of an incoming wave due to the changing mass. Before solving this 
equation of motion exactly, let us first solve for the asymptotic mode functions to gain intuitive understanding. In 
the asymptotic past (t — > — oo), equation (|A7|) is solved by: 

= -«= cxp [-la***] , (A10) 

V2Win 

i.e.: one right- moving or incoming wave with frequency: 

Wi n = (k 2 + A-By . (All) 

In the infinite asymptotic future, the solution necessarily is an appropriately normalised linear superposition of a left- 
and right-moving wave: 

(j)° k ut {t) = a k 2. cxp Huwt] + fa J- cxp [iu) out t] , (A12) 

V 2w ou t V 2Wout 

where: 

iu out = (k 2 + A + BY , (A13) 

and where: 

|K|| 2 -||/3 fe || 2 = l, (A14) 

for a consistent canonical quantisation. In both asymptotic regions we can now immediately calculate the statistical 
propagator from equation (|A8b)) : 

F- m {k,t,t') = 2—cos(oj in (t-t')) (Al5a) 
2w in 

F ont (k,t,t') = ^kla^ + l^cos^it^t'^ + akPie-^^ . (A15b) 

Using equation (|9]), we can calculate the area in phase space the in and out states occupy: 

AL n (t) = 1 (Al6a) 

A^*(t) = 1. (A16b) 

Hence for the entropy we find: 

S£(t) = = Sr\t) . (A17) 
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Figure 17: Causal Propagator with a constant mass. 
Parameters: k/p = 1, m^/p = 1. 




Figure 18: Causal Propagator for a small change in the 
mass. Parameters: k/p = 1, A/p 2 = B/p 2 = 1/2. 




Figure 19: Causal Propagator for a large change in the 
mass. Parameters: k/p — 1, A/ p 2 = B / p 2 = 2. 

We conclude that in both asymptotic regions the entropy i 
mass. 




Figure 20: Statistical Propagator with a constant mass. 
Parameters: k/p = 1, m^/p = 1. 




Figure 21: Statistical Propagator for a small change in 
the mass. Parameters: k/p — 1, A/p 2 = B/p 2 — 1/2. 




Figure 22: Statistical Propagator for a large change in 
the mass. Parameters: k/p = 1, A/p 2 — B/p 2 = 2. 



zero and no entropy has been generated by changing the 
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However, we can do better than study the asymptotic behaviour only. Birrcll and Davies study cosmological particle 
creation in section 3.4 of their book |88( in a simple, conveniently chosen cosmological setting. They consider a scale 
factor as a function of conformal time a{rj) which behaves as: 

a 2 (77) = A + Bt&nh(pri) . (A18) 

This represents an asymptotically static universe with a smooth expansion connecting these two asymptotic regions. 
Indeed, the equation of motion (in conformal time) for the mode functions Birrell and Davies consider coincides 
precisely with (|A7|) . The solution to (|A7|) which behaves as a positive frequency mode in the asymptotic past 
(t — > —00) can be expressed in terms of Gauss' hypcrgcomctric function 2F1: 



y/20J in 



—iui + t — 1 — log{2 cosh(pi)} 
P 



2 F 1 I 1 + 1— , 1— ; 1 - 1— ■ \{1 + tanh(pi)} ) , (A19) 
P P P 2 



such that: 



lim <f%{t) = -== cxp [-tu^t] , (A20) 



where we defined o>i n and w ut in equations (|A11[) and (|A13[) . respectively, and: 



W± = ^(^out ± Win) ■ (A21) 



Alternatively, the modes which reduce to positive frequency modes in the out region are given by: 

■o- w w out 1 

— ,i — ;l + t ;- 

P P P 2 



0rW = ^^exp 

V /w out 



-tto + t — i — log{2 cosh(pt)} 
P 



8^(1 + ,— ,*— ;l + t ^;i{l+tanh(^)}) , (A22) 



such that: 



lim ^ ut (t) = J. exp [-twout*] ■ (A23) 

t^co V2Wout 

We can rewrite the hypergeometric functions using equations (15.3.3) and (15.3.6) of jl!2l | and identify: 

C(*) = a k ct>T\t) + M° k ut *(t) , (A24) 

where 

a k = ( — ^ 3 r( ' 1 - Iu} ^/P) r (- Iu} ° ut /P) (A25a) 

V W in / r (-WJ+//9) r (1 - WJ + /p) 
1 

n _ ( ^out V r (1 - Kjjn/ g) r (zcj out /p) a „, M 

/fc "l^~J r^_/p)r(i + 4W _/ P ) ' (A25b) 

Having the mode functions at our disposal, we can find (the rather cumbersome expressions for) the exact causal and 
statistical propagators. The statistical and causal propagators can however neatly be visualised. Figure fTTl shows the 
causal propagator with a constant mass from equation (|A9a[) for comparison to the changing mass case. In figures [TS] 
and 1 191 we show the exact causal propagator for a relatively small increase of the mass (from m^/p = to m^/p = 1) 
and a larger one (from m^/p = to m^/p = 4) for one particular Fourier mode only (k/p — 1). Figures [20"1 1211 and 
|2"21 show the analogous statistical propagators. 

We can easily relate the statistical propagator to the phase space area by making use of equation ([9]). It will not 
come as a surprise to the reader that we find: 

A fc (t) = l, (A26) 

and hence: 

S k (t) = 0, (A27) 

also for all intermediate times. A final remark is in order. The reader should not confuse |/3fc| 2 calculated from 
equation (|A25b[) with the phase space particle number density or statistical number density (fTTj) . Although the mass 
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Figure 23: This plot shows |/3fc | 2 as a function of the final 
mass mout/p, for fixed k/p — 0.01. The dashed line shows 
the adiabatic regime (m ln _/ p = 0.2) whereas the solid line 
shows the non-adiabatic regime (mi n /p = 0.02). 



is changing, the phase space particle density remains zero but |/3/c| 2 , which in the literature is often referred to as a 
particle number, can change significantly as can be appreciated from figure [231 This is just caused by the fact that 
the in and out vacua differ. Wc plot the behaviour of |/3fc| 2 as a function of m out /p in both the adiabatic regime 
(|/5fc| 2 <C 1) and non-adiabatic regime (|/3fc| 2 ^ !)■ 

This simple example suggests the following: (i) the area in phase space a state occupies is a good quantitative 
measure of the entropy, (ii) the statistical propagator contains all the information required to calculate this phase 
space area, and (iii) a changing mass does not change the entropy for a free scalar field. If we contrast this result with 
the calculations performed in the main body of the paper, it is important to realise that just a changing mass, in the 
absence of interactions, produces no entropy, whereas we have shown that the entropy can change in the interacting 
case. 



Appendix B: Evaluation of the Thermal Propagators 



In this appendix, we will for pedagogical reasons derive the thermal propagators from first principles. The four 
thermal propagators should solve the standard differential equation (fT9f . We start by summarising the conditions the 
thermal propagators have to satisfy in position space: 



iA++(x;x r ) +iA-~(x;x') = iA-+(x; x') + iA+-(x; x') 



iA+ + (x;x')-iA--(x;x') = sgn(t - t')(iA- + (x;x') - iA^ 



~{x;x')) 



iA-+(x;x') = tA+-(x';x) 



*A-+((t- l /3,x);x') 



iA+ (x';x) 



(Bla) 
(Bib) 
(Blc) 
(Bid) 
(Blc) 

Here, the first condition (identical to (|18cl) ) and the second relate the sum and the difference of the time ordered 
and anti-time ordered propagators to the two Wightman functions, respectively. The third condition is jus t identica l 
to (|18dl) . Condition (|Bld|) is the well-known KMS condition or Kubo-Martin-Schwingcr condition, see |l!3l . Ill j j- 
The KMS condition corresponds to periodic boundary conditions in the imaginary time direction due to assuming a 
thermal density matrix operator pth <x cxp[— /3H]. The final equation arises from requiring standard commutation 
relations. Fourier transforming the equations above according to (|38[) yields: 



[ X (t, x), ±(t, x')] = dv {*A-+(.t; x') - iA+-(x; x')}\ t=t , = (x - x') . 



zA++(fc") + *A--(fc^) = iA: 



lA^+ikv) 



iA 
> A 



(*") 
(A") 



-(k-n + *A+~ 
-2t 



k u k^ + ml 



iA-+(k^) 



^A- + (k'")-^A+-(k 



(-ft") 

27rsgn(k°)6{k„k fi +m 2 x ) 



e^iA: 



(B2a) 

(B2b) 

(B2c) 
(B2d) 
(B2e) 
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To obtain the second relation (IB2bl) , we recall: 



f°° 1 
sgn(x) = P / dfc — -( 
J-oo ™k 



(B3) 



where P denotes the Cauchy principal value. Relations (|B2dj) and (|B2cjl trivially yield the two thermal Wightman 
functions: 



tA+-(k») = 2irsgn(k°)5{k^ + m 2 x )n c *(k°) 
iA-+(k») = 2nsgn(k°)5(W + mj)(l + n^(k )) 



(B4a) 
(B4b) 



where n^ q (fc°) is the Bose-Einstein distribution given by (|42|) . In order to solve for the time ordered and anti-time 
ordered propagators, let us make the following general ansatze: 



tA—(fc") = 



kpk^ + m x — le 



knk^ + mz + ie 



+ 6(k^ +m 2 x )f(k°) 
5{k^ + m x )g(k°) , 



(B5a) 
(B5b) 



The functions f(k°) and g(k°) do not depend on k % due to the delta function. We have already chosen the time 
ordered and anti-time ordered pole prescription. This is particularly convenient because, as we will appreciate in a 
moment, this allows us to easily recover the familiar vacuum solutions when T — > 0. We will return to this subtlety 
shortly. Condition (|B2b[) immediately implies: 



f(k°) = g(k°) ; 



where we have made use of the Dirac identity: 



= P iit8{x) 

x + it x 



(B6) 



(B7) 



Because of the time ordering, iA ++ (x;x') = iA ++ (x';x) such that iA ++ (fc M ) = iA ++ (— k^). This consideration 
likewise applies for the anti-time ordered propagator and suggests that the most economic way of writing f(k ) is in 
terms of \k°\. We observe: 



i+n?(*°) = 8 gn(* ) (l +n ^(\k°\)j 
Using the relation above and condition (|B2a[) : 

/(fc ) = 27rn^(|fc°|). 

The thermal propagators are thus given by: 

— i 



iA++(k») 
tA--(fc") 



k^ 



27r£(fc M F 



— = + 2ir5(k^ 

k^k* 1 + mi. + le 



KK q (\k \) 

ml)n<?(\k°\) 



iA+-(k») = 2n5(W+m x )[8(-k )+n x *(\k°\)] 
iA-+(k^) = 2n5(k^ + m x )[e(k ) + n^(\k°\)] . 

If we now let T — > to obtain the familiar vacuum solutions, wc find: 



(B8) 
(B9) 

(BlOa) 

(BlOb) 

(BlOc) 
(BlOd) 



iA++{k^) 



iA x -{k») 



kak^ + m 2 + ie 



iA+-{k^) = 27T(5(fc M ^ + m 2 x )0(-k°) 
iA x + {k») = 2TrS(k fl k" + m 2 x )6(k°) . 



(Blla) 
(Bllb) 

(Bile) 
(Blld) 
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Clearly, writing the thermal propagators in the form (|B10[) above, facilitates obtaining the vacuum solutions easily. 
The reason is that n® q (|fc°|) — > when T — > 0, whereas this statement does not hold for 7i° q (/c°). 

The freedom to choose a different pole prescription such as the advanced or retarded pole prescription, is just an 
equivalent way of writing the thermal propagators, which can easily be verified by making use of the Dirac identity 
(|B7[) . The vacuum and thermal contributions to iA+ + (fc M ) and iA~~(fc M ) separate only so neatly when we use the 
time ordered and anti-time ordered contours to evaluate these propagators, respectively. 



Appendix C: Alternative Method of Renormalising the Self-masses 



In this appendix, we find the correctly renormaliscd self-masses by means of an alternative Fourier space calculation. 
From equations ()34a[) and (|43a[) we immediately deduce: 



zA/++(fc^) = - % — d D (x-x')(iA++(x;x')) 2 e- lk{x - x,) 



(CI) 



ih 2 



d D k' 



1 



1 



(2tt) d k' fc'' 1 + m| - it (fc M - k')(k^ - k'^) + m 2 - it 



-- / dx 



d D l 



(2tt) d (l^ + k^x(l - x) + m 2 - it)' 



where we used Feynman's trick (see e.g. |115| ). performed a Wick rotation and we defined 1^ = k' — xk^ as the new 
Euclideanised integration variable. The integral can now straightforwardly be performed, which yields when m x — > 0: 



iM++(kn 



h 2 (D - 2)(k^ - ie) E 2 ± 2 1 - 2D n a -^ 1 r(^) 
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k^ 
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where in the last line we have expanded around D = 4 as usual and we have again introduced a scale [i to make the 
argument of the logarithm dimensionless. Observe that the (numerical value of the) divergent term coincides with 
equation (f5Tj) as it should be. Note that for iM^~ (k^), we would have to use the other Wick rotation (in order not 
to cross the poles) giving us the desired minus sign difference just as in the position space calculation. Finally, we 
need to perform the k integral in order to derive the self-mass in Fourier space. The relevant integral is: 



32tt 2 
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(C3) 
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In order to make the inverse Fourier integral convergent, we introduce t regulators where appropriate. As an inter- 
mediate result, we present: 



32tt 2 
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I At - it 



I Ail + it 



Clearly, the e regulators in the logarithms and exponents are redundant and can be sent to zero. Using the Dirac rule 
in equation (|B7|) once more, we finally arrive at: 



32tt 2 
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This result is divergent when we let At — > 0. As already discussed in the paper in section UTTl the correct way to deal 
with this is to extract two time derivatives acting on e.g. Z^ + (k, t, t'). We then remain with a perfectly finite result, 
which can be appreciated from equation (1641) . 
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To prove this, let us indeed evaluate the two time derivatives in equation (|64[) . The result is: 
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ic 



The first line contains two elements. The first is obtained directly from evaluating the double time derivative in 
equation (|64[) . The second contribution originates from expanding the term multiplying the delta function in equation 
(|5ip around D = 4, corresponding to the minimal subtraction rcnormalisation scheme. Indeed, the second line of 
equation (|C6[) is identical to (|C5|) as it should be. This shows that the position space and Fourier space calculations 
yield identical results, however the former calculation proves to be superior to the latter as the two extracted time 
derivatives appear naturally in that case. 



Appendix D: Retarded Self- mass 

The retarded self-mass iM^ ren (x; x') can be obtained by means of an independent calculation by making use of 
equation ([52")) : 



iMl iien (x;x') = iM++ n (x;x')-iM+ (x;x') 
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256tt 3 

where as before r = \\x — In Fourier space we find after some partial integrations 
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The last integral is performed in e.g. [11 
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Since the term in square brackets is proportional to (At) 2 as At — > 0, the 6 (At) commutes through one of the (d 2 + k 2 ) 
operators. Evaluating it further yields: 
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If we determine the retarded self-mass directly from equation (|64[) , we find perfect agreement with the result above, 
representing yet another consistency check of ([64]) . 
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Appendix E: The Statistical Propagator in Fourier Space 



Let us now calculate the statistical propagator in Fourier space. Starting point is the equation of motion (|39[) in 
Fourier space. We can straightforwardly write the two Wightman functions as: 



tA+"(fc") = 



where we have made use of the definition of the advanced propagator (|24b|) in Fourier space: 



zAi(fc") = 



k»k»+ml + iMl ien (k») 



Moreover, we made use of: 
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Clearly, we need to evaluate some self-masses in Fourier space. The simplest method of determining e.g. iM^ (fc M ) 
or iM7 (k 11 ) is to use the retarded self-mass in Fourier space and iM^ en (k^) which we already derived before in 
appendix [Cj Using expressions (|T2"j) and (|D4[) we can derive: 
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where fc 2 = (fc ) 2 and where we have used two partial integrations and disposed ourselves of the boundary terms by 
introducing an e regulator where necessary. We can now use: 
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to prepare for another partial integration. Equation (|E4p evaluates to: 
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to evaluate the remaining integrals. The result is: 

h 2 



32tt 2 



log 



-kl + k 2 - isgn(fc°)e 



■2 7E 



(E9) 



From tM^ Icn (k^) = iM£+ n (k fl )-tM£~(k' i ) we can immediately find the Wightman self-masses. We take jM^+J/;'') 
from equation (|C2j) , but we have to make sure we use the same subtraction scheme as in our position space calculation 
in equation ([52"]) . We therefore modify equation (|C2[) slightly to: 
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We thus find: 
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As a check, we consider the following relation that has to be satisfied: 



(E12) 



Of course, iM~f~^ en {k^) is given in equation (|E10[) which also allows us to derive the anti-time ordered self-mass: 
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where iM^ reTl {k^) contains an additional minus sign because of the Wick rotation (see appendix [Cl for details). We 
thus find: 



*M+-{k*) + iM~ + {k») = -j^e(k 2 - k 2 ) . 



(E14) 



This is in perfect agreement with equation (|E1 f p . 

Using equations ([El]) . ()E10j) and (|E11|) we can derive our solutions for the two Wightman functions in Fourier 
space: 
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The limit h — > in the equations above nicely agrees with the vacuum Wightman propagators in equations 
(I43dl) . Hence the statistical propagator in Fourier space reads: 
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